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Abstract 

The HLS Model, equipped with a mechanism providing the breaking of U(3)/SU(3) 
symmetry and an isospin symmetry breaking leading naturally to vector meson mixing, 
has been recently shown to successfully account for e^e^ tt~^it~ cross section and for 
the dipion spectrum in r decay. The present study shows that the full anomalous sector 
of the HLS model can be considered and is validated by the experimental data. Indeed, 
this extended model provides a successful simultaneous fit to the e+e^ —>■ vr+vr^ data 
together with the available data on e+e~ — > n^^, e^e' r/7 and e^e~ — > vr^vr+Tr" cross 
sections. It is shown that the fit of these data sets also predicts an accurate description of the 
rj/rj' ^ 7r+7r^7 decays fully consistent with the reported information on their branching 
fractions and spectra. Finally, one also derives from our global fits products of widths of 
the form T{V h)T(y e+e") and ratios of the form T{V fi)/T{V /2) 
describing decays of vector mesons to several non-leptonic final states. 
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1 Introduction 



A Lagrangian model based on the Hidden Local Symmetry Model 13 has been used 
[[3l in order to account simultaneously for the e+e~ vr+vr" and tt'^tt^Ut data. This 

model has been supplied with a SU(3)/U(3) symmetry breaking mechanism [4','5l, essentially 
relying on the BKY [6J mechanism. In order to apply this model simultaneously to annihilation 
data and to the vrvr spectrum from r decays, an isospin breaking mechanism has been defined, 
which is also a specific implementation of the vector meson mixing. This is needed at least 
in order to account for the to tt+tt^ decay which is an important signal in e+e^ tt+tt" 
data. Moreover, a global fit to e+e^ and r data needs a reasonably general isospin symmetry 
breaking mechanism, active in e~^e~ annihilations and which can be switched off for the r 
decay. 

The peculiar aspect of the mixing model defined in Q is the use of several decay modes 
which serve to constrain the parameters introduced by these breaking schemes. Then, the 
anomalous decay modes of the forrrQ VP^ and the two photon decay widths of the 7r°, i], rj' 
mesons play an important role. 

As such, the model is already overconstrained and has been successfully used |[3l. How- 
ever, it cannot apply to e+e" annihilation channels other than tt+tt^, K^K^ and K^T^ , and 
processes like e+e~ — > 7r°7, e'^e~' 777, e+e~ r]''j or e'^e'^ tt^tt+tt^, remain beyond its 
scope. Indeed, even ifV^ P7 decay widths are an important part of these processes, they do 
not exhaust the physics content of the annihilation processes. 

Therefore, having a model able to describe all these processes (and possibly others, as for 
instance e+e^ — > n^cu) is certainly valuable. This is already important from a physics point of 
view but, more technically, this may allow to check the relative consistency of various data sets 
within the same framework and, therefore, the systematics affecting the data sets. 

An extension of the HLS model is possible by introducing the anomalous Lagrangian pieces 
identified long ago ^ 13 and relying also on the above mentioned breaking schemes reminded 
or defined in |l3l. This could allow interesting improvements in estimating the hadronic con- 
tribution to the muon anomalous magnetic moment a^. Indeed, as basically several crucial 
physics parameters are common to all low energy annihilations, considering all annihilation 
processes altogether may turn out to increase the statistics and, then, improve uncertainties. 
In this way, one may expect a better determination of e+e^ tt+tt^ using all other anni- 
hilation processes. Conversely, the huge experimental effort invested in order to get accurate 
e+e^ TT+TT^ data may help in better estimating contributions to from other processes, 
e+e^ Ti^TT'^n' for instance. However, the systematics present in all these processes may 
limit the expected improvements. 

The purpose of the present study is to construct such a model and test its ability to account 
for the above mentioned processes beside e+e^ tt+tt^. This study will be split up into 
two parts for ease of reading. In the following, we will built the model and apply it to several 
annihilation processes and decay width information. An accompanying paper is devoted to 
analyzing the dipion spectra in r decay and some non-perturbative contributions to a^. 

' Throughout this paper, we denote by P and V, resp. the basic pseudoscalar and vector meson nonets. The 
explicit form of their matrix representation can be found in several papers lH] |2] [H [5J and is reminded in [3J. The 
electromagnetic field is denoted by A. All definitions and notations in the present paper closely follow those in 
0. 
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We leave aside from the present study the e+e^ K^K^ / K'^K annihilation channels 
which are known to raise a problem thoroughly examined in [8] and may call for a solution 
beyond the scope of effective Lagrangians flU. Indeed, the ratio of the partial widths of (p 
decays to a charged or neutral kaon pair looks inconsistent with expectations. This question 
raises an issue which deserves a specific study at the level of cross sections which we have only 
started. 

The paper is organized as follows. Section [21 reminds the ingredients to be used in order to 
construct the Lagrangian of our Extended Model : these can be found in [|71|2l. Without further 
constraints, this leads to introducing three more free parameters. In Section[3l we require these 
anomalous Lagrangians to recover identically amplitudes at the chiral limit as given by the 
standard Wess-Zumino-Witten (WZW) Lagrangians [fTOl[TT| ; this leads to constraining some 
of the parameters above and, then, to reduce the additional parameter freedom of the model. 
In Sections |4] and |5l we briefly remind the isospin breaking model of [|3l, mostly in order to 
simplify notations. Sectionals devoted to our treatment of correlated systematic errors which 
introduce uncertainties on the global scale of spectra ; we also emphasize the way to check 
the scale information provided with the various data sets. A method to account for possible 
data discrepancies is also given. In Section |71 the Extended Model is shown to keep intact the 
description of the e+e^ tt+tt^ data examined in Q, together with the unavoidable set of 17 
radiative and leptonic decays. 

In Section [U the real use of the Extended Model begins with fitting the e+e" — > 7r°7 and 
e^e^ r]'j data simultaneously with our usual set of e+e^ tt+tt^ data fS^. While introduc- 
ing these data sets, one has to remove all Vn^'j and Vr]'j decay pieces of information from the 
set of decay partial widths used together with annihilations, in order to avoid redundancies. At 
this step, the new Lagrangian pieces are shown to correctly account for the data. In Section |9l 
one performs a fit of the e+e^ — ^ tt+tt" data collected by the KLOE Collaboration using the 
Initial State Radiation Method (ISR). This is seemingly the first published fit to the KLOE data 
[fT2l[T3l . The next step, including the e^e^ tt^tt+tt^ annihilation data in our fit procedure, is 
the subject of Section[10l where a discussion on the various data samples is done. In Section [TT] 
numerical values for to and meson parameters are given and the exact strength of the APPP 
and VPPP couplings is emphasized. 

The aim of Section [l2l is to show that the Extended Model, with parameters fixed by fitting 
to the cross section listed above, allows to accurately predict the properties (spectrum lineshape, 
partial width) of the decay channels rf/r]' ^ 7r+7r^7, both dominated by the box anomalies of 
QCD. In Section [131 we focus on vector meson partial widths and discuss the question of 
getting their absolute values from e+e^ data only, where models are sensitive - for each vector 
meson - to only the product of its leptonic decay width and of its partial width to the final state 
of the e+e^ annihilation. Finally, a few concluding remarks are collected in the Conclusion to 
the present part. Some reminders or results are given in Appendices, in order to keep the main 
text as clear as possible. 
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2 The Extended HLS Lagrangian Model 



The most general Lagrangian of the HLS model El 13 involves anomalous and non- 
anomalous sectors. It can be written : 

C- = ChLS + ^YM + ^WZW + ^FKTUY (1) 

The piece denoted Chls can be written Ca + aCy, where a is the specific HLS parameter, 
expected equal to 2 if vector dominance is fully satisfied. This Lagrangian piece allows a direct 
coupling of photons to pseudoscalar mesons to survive with a magnitude proportional to a — 2. 
This piece, which accounts for the non-anomalous sector, has been given in expanded form in 
|4J. A brief survey of the HLS model can be found in Appendix A of [3 J and is enough for 
the purpose of the present paper. A recent and comprehensive review can be found in [2], for 
instance. 

The parts of the HLS Lagrangian Chls specific of the e+e" physics (Cvmd) and r decays 
(Ct) have also been given in [|3|| ; they are reproduced with almost no comment in the Appendix 
(see Section |A]) in order to avoid too frequent references to an external paper. 

Cym is the Yang-Mills piece which accounts for the vector mesons and exhibits the non- 
abelian group structure of the vector fields (see Eq. (D.l) in jSj). 



Cwzw and Cfktuy are the anomalous Wess-Zumino-Witten [fTOl [TTIl and FKTUY 
lITllH Lagrangian pieces which account for parity violating decays. Very briefly, the anomalous 
Lagrangian can be formally written : 



^anomalous Cy/ ZW ~l~ ^ ^ Q Hi 



(2) 



where the Cj's are arbitrary constants weighting the additional FKTUY anomalous Lagrangians. 
Limiting oneself to the photon and vector meson couplings, Canomaious can be cast in the 
formO : 

^anomalous = ^VVP + ^AVP + ^AAP + ^VPPP + ^APPP (3) 

where A denotes the electromagnetic field. The Lagrangian pieces occuring in Eq. ([3]) areH 
Ul : 



Cvvp = 

^AVP = 

Caap = 

^VPPP 
^APPP 



Kg' 

Ncge 
Nce^ 



(C3 - c,)e>''"'f'd,A,Ti[{d^Vf,, Q}P] 
(l-c,)e'^'^-''d,A,d^ApTT[Q^P] 



(4) 



Svr^/^ 4 



-For clarity, the new constant parameters are denoted exactly as they are defined in [ID. 
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where we have extracted the universal vector meson coupling constant g from the definition 
of the vector field matrix (see [3J for our definitions), Q denotes the quark charge matrix and 
{daVi3,Q} is the anticommutator of its arguments. Nc is the number of quark colors, 3 in our 
case. The original WZW terms here can be easily identified in Eqs dH) and one also observes 
that these Lagrangians depend on only 3 (FKTUY) parameter combinations : Ci — C2, C3 and 
C4. We may refer to the first three Lagrangians in Eqs. (U) as Triangle Anomaly and the last 
two as Box Anomaly Lagrangians. 

If full Vector Meson Dominance were fulfilled, one would have : 

C3 = C4 = 1 and ci - C2 = . (5) 

Indeed, one can easily check that these conditions turn out to cancel the 3 Lagrangian pieces in 
Eqs. (HI) involving the photon field (A). In this case, photons connect to hadrons only through 
the V — A transitions allowed by the non-anomalous HLS Lagrangian (see Appendix lAl). 

However, relying on a limited subset of data, it has been argued in |[71 l2l that some experi- 
mental data may instead favor : 

C3 = C4 = Ci - C2 = 1 (6) 

Previous studies especially devoted to VP^ and P77 partial widths f5] performed assuming 
Eqs. ^ were succesfuU enough to conclude that these conditions were well accepted by the 
data. Our work lfT4l on a simultaneous account oi rf/r]' — > 77 and rj/r]' ^ ix^-n^'y led to the 
same conclusion. This gave a strong support to the relevance of the box anomaly phenomenon 
in pseudoscalar meson decays. 

Therefore, the present status is that full VMD is well accepted by all existing data only in 
the Triangle Anomaly sector, while the rf/r]' — > 7r+7r^7 decay modes prevent extending full 
VMD to the Box Anomaly sector. 

However, the statement C3 = C4 = ci — C2 = 1 cannot be considered as firmly established 
without being motivated by a detailed study of all relevant data. If any, such a conclusion should 
only follow from a study involving the annihilation processes e+e^ 7r°7, 777, rf'^y, tt+tt^tt^ 
in their full intricacy, beside the properties of the decay modes ///r^' ^ 77 and ri/r]' ^ 7r+7r^7. 
As anomalous transitions also involve non-anomalous pieces, the data set to examine should 
certainly include the e+e^ tt+tt^ annihilation process in order to define precisely the other 
model parameters. 

Whether such a program is realistic and can be practically worked out is the subject of the 
present paper. It implies to go much beyond what has been done recently in [3] and formerly 
in studies devoted to radiative decays, our [ISllMl for instance. 

3 U(3)/SU(3) Breaking And The WZW Conditions 

We will not rediscuss here the U(3)/SU(3) breaking conditions in the non-anomalous sector 
which have been extensively discussed in our flU |3l following the idea of Bando, Kugo and 
Yamawaki [61, referred to as BKY mechanism. A detailed summary of the U(3)/SU(3) (BKY) 
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breaking scheme is given in Appendix C of ['S^ where the corresponding breaking for the Yang- 
Mills piece has been also considered (see Appendix D herein). Let us only remind that the 
pseudoscalar bare field matrix P (= Pg + Pq) undergoes a renormalization to P' : 

P'^P^ + xP^ = X'f{Ps + Po)X'f 

(7) 

^ Xa = Diag(l, 1,za) 

where x ^ 1 carries the (U(3)) nonet symmetry breaking!^. The other breaking parameter has 
a well defined expression in terms of the kaon and pion decay constants : za = [/x/Ztt]^- 
Correspondingly, the renormalized vector field matrix V is given in terms of the bare one V 
ll3Jby : 



V 



^1/2.. yl/2 



(8) 

Xt = Diag(l, 1, zt) 

Considering a HLS Lagrangian as general as allowed from its derivation [|7l |2l| leads to 
introduce a priori three more basic parameters (ci — C2, C3 and C4) beside the standard g, a (and 
e) occuring already in the non-anomalous Lagrangian. 

However, one may impose a natural constraint to this Extended HLS Model by requesting 
that it provides amplitudes identical to those derived from the WZW Lagrangian alone at the 
chiral point. These amplitudes may only be modified by U(3)/SU(3) breaking effects affecting 
the pseudoscalar sector [l5l[3]|. Instead, additional SU(3) breaking effects which may affect the 
trace terms in Eqs. @i, as well as the dependence upon the c^'s, should drop out at the chiral 
point. 

We perform the SU(3) breaking of each Lagrangian in Eqs. (HI) using a prescription inspired 
from Bramon, Grau and Pancheri ifTTlfTSl . This motivates the following breaking scheme : 



Cvvp '■ 


Tild^V^d^VpP] =^ Ti[d^V,Xwdo.VpP] 


Cavp '■ 


d^A,TT[{d^Vp,Q}P]-- 


^d,A,TT[{d^Vp,QXu}P] 


* C-AAP '■ 


d,A,d^ApTi[Q'P] = 


> d,A,d^ApTT[XnQ'P] 


C-VPPP '■ 


TT[V,d,Pd^PdpP] = 


> TT[XKV,d,Pd^PdpP] 


. Cappp '■ 


A^TT[Qd,PdaPd(,P] -- 


=^ Ti[XLQd,PdaPdf3P] 



(9) 



This might introduce as many as 5 more (breaking) parameters (zw, zu-, zji, zk, zl) as 
the five matrices introduced here are supposed to carry the same form than Xa or Xt reminded 
above. We now show that, actually, it is not the case. 



3.1 The Triangle Anomaly Sector 

Using the Lagrangians in Eqs. (HI) broken as shown in Eqs. ^ just above and the 'j — V 
transition amplitudes as given by the non-anomalous Lagrangian, one can derive the ampli- 

^See ifTSl for the relation between this parameter and determinant terms [16] explicitly breaking the Ua(1) 
symmetry. 



5 



tudes : 



One gets first : 



T 



X 



77) e^'^'-^qUleyp 



Tin' ^ 77) 



a 



1 + 2(C3 - C4)] 



7r°, ?7, r/') 



(10) 



(11) 



which implies C3 = C4 in order to recover the usual WZW term. Therefore, the Cavp La- 
grangian vanishes. Using this condition, one similarly derives : 



T{m 77) = -X 



T(r/8 77) 




- C3){zwZt - Zr) ZwZt + 5za 



3za 

-2(1 - C3){zwZt - Zr 
"iZA 



+ 



+ 



'iZA 
hZA — 2z\yZ^ 

3za 



(12) 



for the singlet and octet parts of the rj and rj' mesons. In order to recover the usual WZW 
expressions (see for instance |5ill4l or more recently [3J), one clearly needs to request zwz^ = 
1, which was first phenomenologically found as a numerical constraint arising from fits fS^. If, 
additionally, one requests zr = 1 (which may look quite natural), then the standard WZW 
amplitudes for r^o/s 77 recovered without requiring any further constraint on C3. We are 
thus led to choose as constraints on the parameters : 
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ZwZj' 



zr 



C3 = C4 



(13) 



3.2 The Box Anomaly Sector 

As measurements related with box anomalies only involve couplings of the form Xn'^n^'j 
with X = 7r°, 1], r]', we only focus on this sector. Let us list the Lagrangian pieces relevant for 
this purpose. 

• A part only of the VVP Lagrangian plays a role in these couplings : 



C 



VVP 



V8 



+ 



.2^3 VQ 



+ 



ZiyZrp 

za 



V8 _^ xrjo 
73 VG 



(14) 



Nc9^C3 
in terms of ideal fields. 

Actually, the ^VP couplings relevant for our purpose can be derived from the following 
effective piece constructed from Eq. (fT4l) above on the one hand, and from the — V 
transitions of the non-anomalous Lagrangian [[3l on the other hand (see Appendix |A] 
below) : 



AVP 



F = 



t7T 



-V8 



C' = 



NcQ ec3 
— vr 



1 



X 



Za 



X 



;V8 



(15) 
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For the purpose of constructing the Xn^n^j coupling at the chiral point, this piece fully 
replaces Eq. (fT4l) . It will be used together with the transition — > tt+tt^ given by the 
non-anomalous Lagrangiar0. 



The VPPP piece of relevance is : 



CvPPP 

D = 



Ncg 



Pi 



12 

C1-C2- C3) 



duVs 



X 



n/6 
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The relevant APPP piece writes 

1 



^APPP 

E = 



4 '12 ~vr 



dnTX d, 



1 - ^(Cl - C2 + C4) 



(16) 



(17) 



With these pieces at hand, one can compute the amplitudes at the chiral point. One first 



gets : 



T(7r° — >• TT^TT 7) 



ieNr 



1 + ^(C3-C4) 



(18) 



which coincides with the amplitude expected from the WZW Lagrangian alone if C3 = C4 
Assuming this condition, one can easily derive at the chiral point : 



T{t]s TT + TT 7) 



ieNr 



T{r]o TT+TT 7) 



12^V3 



X 



(19) 



which proves that the WZW usual (U(3)/SU(3) broken) amplitudes [fT4ll at the chiral point are 
recovered. One may notice that the breaking matrices and Xi play no role in the Xtt+tt^7 
sector. 

Therefore, assuming the conditions summarized in Eq. (fT3l) leaves us with 2 more physics 
parameters unconstrained, compared with the previous version of our model : C3 and ci — C2. 
Using the information previously defined, the Extended Lagrangian we propose is : 

C = Ca{za) + aCvizv) + Cym{zt) + 

+ Cvvpizw) + ^aap{zr = 1) + Cvppp{zk = 1) + ^APPpizL = 1) 



(20) 



C4 = C3 



. ZwZt 



Anticipating on what is reminded in the next Section and relying on (3], this is justified by the fact that ideal 
vector fields and physical vector fields coincide at s = 0. Then, at s = 0, the couplings of the physical oj and (j> 
mesons identically vanish, s being the square of the momentum carried by the vector meson. 
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where each SU(3) breaking parameter is exhibited with its value (when fixed). Out of these, 
the breaking parameters to be determined from data are za-, zy-, zt, x, keeping in mind that 
za = [fx/ fn]"^ should always be satisfiecj^. 

Such a price for the extension of the HLS model we propose looks acceptable, taking into 
account that we plan to describe all available anomalous processes beside e+e^ tt+tt^ and 
t"^ — > TT^'^Tr^z/,- already accounted for within this model. This covers the e+e^ 7r°7, 777, tt+tt^tt' 
cross sections and the box anomalous processes 7]/r]' 7r+7r^7. A priori e+e~ — > r]''j data 
will fall into this extended scope when they will become availablel^. 



4 Isospin Breaking And The (p , tJ, 0) Mixing 



We adopt here the {p ,uj, 0) mixing developed in [3 1 . We only remind here a few properties 
for completeness and to simplify the notations. Let us define : 



Ms)) 



e2{s)=glj,AU+{s) + Uo{s)) 



(21) 



9',..T^'is) 



where n'(s) denotes the tt+tt^ amputatec^ loop function, n_|_(s) and no(s) the charged and 
neutral amputated kaon loops. The coupling constants occuring in Eqs. (|2TI) fulfill gp^^^, = 
'2zAgpKK = 05'/2, in terms of the basic parameters of the model and of the SU(3) symmetry 
breaking parameter za = [//^//tt]^- These loop functions are analytic functions each real on 
the real s-axis below the corresponding threshold. At the limit of equal charged and neutral 
kaon masses, ei(s) vanishes ; on the other hand, both ei(s) and £2(5) have small magnitudes [3j 
in the whole s-region we are interested in (from the 2-pion threshold to the mass). Instead, 
n^7r(s) is basically the p self mass and is thus known to be significant in this s-region. The 
functions in Eqs. (|2TI) are given by dispersion relations and contain each a polynomial in s 
chosen [3 1 of second degree and vanishing at the origin. These functions mainly serve to define 
three complex quantitieq^ : 



( a{s) \ 

m 

\ lis) ) 



n^7r(s) - 62(3) 
pei{s) 



[1 - zv)m? + n^^(s) - p^e2{s) 



pe2[s} 



V 



(l-Zy)m2 + (l-/i2)e2(s) 



(22) 



^It could even be directly fixed to its experimental value [Jk //tt]^ — 1.495 ± 0.031 without degrading the fits, 
process like e+e^ ujtt^ is of this kind, however, one has first to carefully study the effects of scalars 
before any use of the corresponding data. 
^i.e. with unit coupling constants. 

^ The functions a{s), (3{s), 7(5) defined here may be used all along this paper without exhibiting their func- 
tional (s) dependence. This notation makes easier reading the formulae given in Section 8 of |i3J which are used 
in the present study. In order to avoid confusion, the fine structure constant is always denoted aem- 
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where zy is a parameter involved in the SU(3) symmetry breaking of the Cy part of the non- 
anomalous HLS Lagrangian [l4l[5l[3l. We have defined fi = zy\pl, and vn? = ag'^f^ is the 
Higgs-Kibble p meson mass squared generated by the spontaneous breakdown of the HLS 
model. The complex quantities a{s), P{s), 7(s) are, in some sense, "angles" which describe 
the (p , UJ, 0) mixing. At lowest (first) order in ei(s) and £2(5), the mixing scheme is given by : 



UJ 



f p° - auj + 13(f) \ 

UJ + ap^ + 'jcj) 
\ 4>- f3p^ -juj j 



(23) 



These relations exhibit the connexion between the "ideal" vector fields (carrying definite 
isospin properties), which are entries of the vector field matrix, and the "physical" fields (each 
a mixture of isospin and 1 fields) which enter the physics processes. One may consider that 
a{s) is the p — uj (complex) mixing "angle" and that /5(s) and 7(3) are, resp. the p — 4> and 
UJ — (j) mixing "angles". In our model, this mixing is s-dependent ; for instance, the uj — cj) 
mixing angle 7(5) has not the same value at the uj mass and at the (j) mass. Moreover, among 
these angles, 7(5) is the single one to be practically real up to the (p mass region. Finally, as a 
consequence of a(0) = /9(0) = 7(0) = 0, physical (s-dependent) and ideal fields coincide at 
s = 0. 



5 The p Propagators And The 7 — y Transition AmpHtudes 

The inverse propagators for the charged and neutral p can be written : 

r DpQ{s) = s-m^ -Yipp{s) 

(24) 

[ Dp± (s) = s - m2 - 5m^ - Wpp{s) 

allowing for a possible p^ — p° mass difference. Neglecting the effects of mass differences be- 
tween charged and neutral pseudoscalar mesons in the propagators, the same self-mass occurs, 
approximated by [|3l : 

Iipp{s) = Y\'pp{s) = n^^(s) + €2(5) (25) 

The notations defined in the previous Section allow one to express the 7 — F transition 
amplitudes in a more readable way than in [3J (see Section 8 there). The 7 — V transition 
amplitudes from photon to physical vector fields can be written as eFy^{s) with : 

Fy-y{s) = fy^ - ny^(s) , (26) 

the constant term can be read off the {V ■ A) terms in the non-anomalous HLS Lagrangian and 
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the loop correction depends on the pion and kaon loops. The constant terms become 

Uy = agfn 

and the s-dependent loop terms are (31 : 

^ QpTTTV ^ QpTTTV 



1 fl ■ 

a H — 7 

3 3 . 



'1 _ ^^ ^Injs) ^ _ a _ ^^ ^ijs) + €2(3) ^ ^ 62(3) -ei(g) 

2 QpTTTT 2 QpTTTV Q pTTTT 



n 



n 



9pn; 

mjs) 



(1 - 



(^A - - - &)/^ 



QpTTTV 

ei(-s) + ^2(3) 

QpTVTT 



QpTVTT 

£2(3) -ei(3) 

QpTVTT 



(27) 



(28) 



with /i = 2;y-\/2 and b = a{zv — l)/6. n^7r('S) (see Eq. (l25l) ) and HZ-n-i^) carry different 
subtraction polynomials [3J. In the fit procedure described below, as in [[3]|, their (second 
degree) subtraction polynomials are chosen independently and fit from data. Numerically, one 
finds no significant correlation among these 2 polynomials. 



6 Combining Statistical and Systematic Uncertainties 

For any of the data sets we use, there are reported statistical and systematic errors. One 
way to proceed is to add them in quadrature and define correspondingly a to be minimized. 
If the errors are large enough, there is no real need to go beyond this simple treatment. 

However, in samples where statistics is large, systematics should be handled a little bit 
more carefully. Systematic errors can be split up into two different kinds : uncorrected and 
correlated uncertainties. It is quite traditional to combine uncorrected systematic errors and 
statistical errors in quadrature and we follow this rule as, moreover, several data sets are pro- 
vided with this combination already performed. When statistical errors are mentioned, this 
combination should be understood, unless explicitly stated. 

It remains to handle the correlated systematic errors. As a first statement, one may interpret 
these as reflecting a global scale uncertainty which affects the data set considered. This was 
already done in our previous analysis [[3l. Let us rephrase it with slightly more details, in order 
to explain clearly the method which underlies the present study. 



6.1 Scale Uncertainties, A Reminder 

The way statistical errors and scale uncertainties combine can be treated rigorousl>|^. Let us 
assume one has a data set m : {m,, i = 1,- ■ - n} and, correspondingly, a model M : {Mi, i = 

'We gratefully ackowledge P. Astier, LPNHE Paris 6/7, for several dicussions on this subject. 
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1, • • ■ n}. Let us assume given its (symmetric) statistical error covariance matrix V which needs 
not be diagonal. Let us finally assume that a normalization scale uncertainty A affects the data ; 
A is supposed to have as most probable value and standard deviation a. 
In this case, the conventional is : 



= [m-M- AXf [m - M - AX] + '— (29) 



A^ 



a 

where A is traditionally the vector of the model values M and the other notations are obvious. 
For definiteness, in our study, we preferred using A = —m in order to avoid having a covariance 
matrix depending on fit parameter^ 

In this approach, A is nothing but an additional constrained fit parameter. One can solve 
this equation at minimum for A using dx^ /d\ = and substitute the expression for A in Eq. 
(|29l). This leads to : 

x'^ = [m-Mf W-\a^) [m - M] (30) 

where : 

W-\a^) = [v + a^AA-]-' = - rT^I^^^V^^^"^)^^"^)' ^''^ 

which corresponds to Eqs (46) and (47) in ['S^. If the expected value of A were some Aq ^ 0, 
one has to replace in the second term of Eq.(|29l) A^ by (A — Aq)^ and m — M would become 
m — M — AXq in Eq.(|30l). For practical use, if a scale uncertainty has been identified, data are 
generally corrected for this and then Aq = is justified. 

Eq. (|30l) illustrates that a correlated scale error is algebraically related with the model : 

A^^^^"[--f (32) 
ATV-^A+ — 

If the model depends linearly on parameters to be determined, the substitution has certainly 
to be performed in order to avoid the error covariance matrix of the (fit) parameters having a 
zero eigenvalue. If the dependence is non-linear, avoiding solving for A only increases errors 
and produces (spurious) correlations. 

Dealing with one (or several) data sample(s), the value of A following from minimizing 
Eq. (|29l ) can be confronted with the expected correlated systematic uncertainty. In practice, if 
the mean correlated systematic error has its correct value, one expects \\^''^\/a small enough. 
In this case, having checked that A = is reasonable as expected, one can impose from start 
A = 0, i.e. minimize Eq. (|30l) ; one should then check that the fit probability and the parameter 
(hidden inside M) central values are nearly unchanged while the magnitude of their errors 
decreases. 

In lO, we performed slightly differently : the scale uncertainty was considered as a random 
distribution 5\ of zero mean and of standard deviation a which affects the measurements. One 
can check that this approach leads to the same conclusion which is summarized by Eq. (|30l) . 
The present approach only clarifies that the final fit should be performed with A = 0. 



We checked with several of the fit configurations described below the difference between the two possible 
choices A — M and A — —m. We did not observe differences beyond the 0.3 ct level for the fit parameter values. 
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6.2 Checking And Dealing With A Missing Variance 



However, as well known, identifying and estimating systematic uncertainties can be a del- 
icate matter. Some source of systematics could have been missed or underestimated. When 
dealing with only one data set, this could well be quite transparent, as its effects could be 
absorbed by the other parameter values. However, while merging different data sets with dif- 
ferents systematics, the pattern can be quite different and could result in poor fit probabilities. 
Therefore, it is useful to be in position of identifying a possible missing variance affecting some 
scale. 

Let us assume the correct systematic variance of A be cr^ + rf instead of cr^. In this case the 
correct formula is Eq. (|29l ) with o"^ ^ a'^ + 7f. The question is now how to detect that a piece 
represented by rf could have been missed. 

One can first check that : 

= [m - M - A(Ai + V-^ [m-M- A{\^ + A2)] + 4 + 4 (33) 

allows to recover the right result, provided one treats Ai and A2 as independent variables of 
zero mean and of respective variance cr^ and rf. One can, moreover, check that : 

= \m-M - A\2f W-^ (a^) [m-M - AA2] + (34) 

by solving, just as before, for only the identified systematics represented by the couple (Ai, a^). 
This equation gives us a handle to account for a missing (part) of the variance which could be 
revealed by using a large ensemble of data sets, each with its own systematics. 

A way to check for a possible missing variance, is to compare a fitted A2 to the "identified" 
variance cr^. A tentative assumption for rf could be to state rf c:^ a'^. Then, within a numerical 
fit procedure, one can check the magnitude of ja. If |A2**|/o" is small enough (< 1^2) 
a missing variance is certainly negligible compared with the identified one. Otherwise, the fit 
value of the scale indicates how much the data should be rescaled in order to match all other 
data sets and the model under test. 

Of course, one thus makes an implicit statement : if the data sets which exhibit a significant 
missing variance represent a small minority, one may consider that this validates both a correc- 
tion for missing variance and the model. If, instead, the data sets exhibiting a missing variance 
represent a majority, the model is certainly invalidated. Intermediate situations, if any, would 
be uncomfortable ; this might indicate some unaccounted for physics. Anyway, the set of data 
samples we consider does not face us with the latter configuration. 

Finally, the we shall deal with is a sum of partial of the kind shown by Eq. (|34l) . A 
first run of the fit procedure helps to identify which of the A2 (a being the data set index) can 
be safely dropped out. 

Concerning those which have still to be considered, one can keep the fit parameter as intro- 
duced. One should however keep in mind that this certainly enlarges the variance. 

As a final remark, one should mention a rigorous way to lessen the variance. If some A2 is 
considered significant, one should account for it, as the corresponding data set has certainly not 
been corrected for this source of uncertainties. However, a numerical minimization procedure 
like MINUIT provides an estimate of A2 - we name it A2** - and an estimate of its uncertainty 
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- we name it a^j^g - which can be accurately known running a code like MINOS. Then, the 
corrected (partial) to be minimized can be rewritten : 



X 



m-M - *A 



m-M - Xi'^A 



(35) 



can be estimated from Eq. (1321) by changing to rj'^ = cTmiss- '^^^ removes one free 
parameter, one almost certainly lessens the variance. This last expression depends on the other 
parameters under fit. A quite acceptable solution is instead to use the numerical minimizer 
outpu10 for both and cr^j^^ , which makes the convergence easier. 



7 The First Step In Modeling 

In order to compare the pion form factor in e+e^ annihilations and in r decays, one has 
to account for isospin breaking effects which differ in both processes. The purpose of [3] 
was to show that the mixing between the to and (p mesons is responsible for most of the 
reported difference. We succeeded in determining the mixing model sketched above (mostly 
the "angles" a{s), P{s), 7(s)) using the information provided by the anomalous decays of type 
VP'j (and P77) and the information carried by the isospin violating decay^ cu/^ — ^ rni. In 
order to constrain more efficiently the parameter set, the partial widths V e+e~ (V = cu, (p) 
were also included. These pieces of information, as well as the co/cp mass and width, were 
fixed at their accepted values [fT9ll . We plan to examine the behavior of the Extended Model 
presented above while introducing more and more information to account for. For this purpose, 
we first focus on e+e^ data, and as a first step on the e+e~ tt+tt^ data mostly collected at 
Novosibirsk and already examined in [3J : 

• the former data sets collected by the OLYA and CMD collaborations [l20l, and by DM1 
at ACO fi2ril ; these were (and are still) together referred to as "old timelike data", 

• the data sets more recently collected by the CMD-2 [|22l|23l|2ll25l and SND [26l Col- 
laborations, referred to globally as "new timelike data" 

• all the partial widths for the decay processes of type VFj, P77, u/cp — e+e^ and 
— > TT+TT^ at their updated recommended values [fT9l . These represent 17 pieces of 
information. 

We gave up including the decay width rj' p7 (actually i]' 7f~^TT~^) as, in addition to 
the dominant triangle anomaly contribution, there is some (small) contamination by the box 
anomaly discussed already above. The decay modes rf/r]' ^ 7r+7r^7 are revisited at the end of 
this study. 

The systematic errors on the e^e" tt+tt^ cross sections just quoted are treated exactly 
as explained in Section 11.2 of [3], or as reminded above in Section |6l i.e. by summing in 
quadrature the statistical errors and the uncorrelated part of the systematic errors on the one 
hand, and, on the other hand, by accounting for the correlated systematic error through a global 

"in our case, the non-linear character of the model avoids having a singular parameter error covariance matrix 
while keeping the scale among the parameters to be fit. 

'^Actually, the w — > tt+tt^ mode has not to be included as it is already part of the e+e tt+tt^ spectrum. 
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scale to be fit, at least for cross-check, as also explained in Section [6l The corresponding 
standard deviations to be introduced in the are 0.4% for the "new timelike data" and 1.0% 
for the "old timelike data" ETl . The fitted scale values are expected negligible 

As our data samples contain all data used in order to get the recommended mass and width 
for the uj and mesons [T9l, we start here by leaving free the uj mass and width ; this could 
have easily been avoided at this step of our study by using [3] the world average values. 

In addition, the present work uses the code delivered by F. Jegerlehner in order to calculate 
Ciem{s) wWch Can be downloaded from [|28ll ; this code, partly documented in ||29l , has been 
constructed in order to improve the estimates of the muon g—2 and of aem(s) [l30l[3Tl[32l . From 
there we get the photon (hadronic+leptonic) vacuum polarization (VP) building block, which 
is used by multiplying our model form factors, generically denoted F{s), by the corresponding 
factor : 

ns) =^ , , , F{s) (36) 

The package alphaQED.uu [|28l also provides an estimate of the uncertainty on Aaem(s)- 
In the extended model of this paper, the pion form factor has the same expression as in []3l : 



S] 



(37) 



where Dp{s) is defined by the first Eq. (|24|) and by Eq. (l25l) . while the other propagators are 
standard fixed width Breit-Wigner formulae. With our new notation set, the coupling constants 
simply write : 

ag ag ag 

gpTTTT = — , g^^-K-K = — ^a(s) , g^nn = yPl-s) (38) 

On the other hand, the VP'j coupling constants entering the corresponding partial widths 
in the extended model are essen tiall>@ given by those in Eqs (E.1-E.4) of [SJ multiplied each 
by the new fit parameter C3. Instead, the couplings constants of type P77 given by Eqs (E.5) in 
[m are left unchanged as well as the leptonic decay widths of the vector mesons. 



For definiteness, we have first performed the fit without including the photon VP. The 
most interesting results are reported in the first data column of Table [U The results obtained 
when introducing the photon are given in the second data column in the same Table. One 
may already conclude that both descriptions provide a quite good account of the data sets 
considered. The first data column in Table [21 displays the fit value of the parameters having 
the most intuitive meaning while fitting with the photon VP. The values found for the fit scale 
factors exhibit a nice correspondence with the expectations reported in the experimental papers. 
One should also note that the fit value of the newly introduced fit parameter C3 is statistically 
consistent with 1. This is in good correspondence with the fits presented in [[31. 

We do not discuss any further the intrinsic fit quality and fit parameter values as this case 
becomes interesting only compared to what happens when using additional data sets. We also 
do not show plots illustrating this fit quality : they are visually indistinguishable from Figure 
2 in [3]. One may, however, remark that the probabilities are more favorable now ; this should 

'^Negligible correction terms are outlined in Section[T3] 
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Data Set 
ft (data + conditions) 


Without VP 
NSK 


NSK 


Vith Vacuum ] 


polarisation (\ 
++ KLOE 


^P) 

+++ 7r°7r+7r" 


Decays 


7.78/(17) 


7.77/(17) 


14.31/(9) 


14.60/(9) 


14.70/(9) 


New 
Timelike (127+1) 


114.83 


114.09 


114.12 


127.15 


127.75 


Old 

Timelike (82+1) 


54.80 


53.86 


50.84 


49.32 


49.32 


KLOE (60+5) 








108.39 


108.21 


TT^ (86) 


— 


— 


61.38 


62.07 


65.66 


777(182) 






128.55 


129.73 


135.20 


xVdof 
Probability 


177.41/210 
95.0 % 


175.72/210 
95.9 % 


369.50/468 
99.9% 


491.27/528 
87.2 % 


637.90/653 
65.6 % 



Table 1: The first data column displays the partial information while working with the 
largest set of partial widths and all available e+e^ tt+tt^ data sets (except those of KLOE) 
and not taking into account the photon VP. The second data column displays the corresponding 
information while, instead, introducing the photon VP effects. In the third data column, the 
data concerning e"^e^ 7r°7 and e~^e~ 777 are considered together with those on e+e^ 
TT+TT^. In the following data column, the KLOE data are included and, correspondingly, the last 
data column gives the fit information while adding also the e+e^ tt^tt+tt^ data (see text). 
Boldface numbers in the first data line display the number of (independent) partial widths which 
are included in the full fitted data sets. 

be attributed marginally to using a different photon hadronic vacuum polarisation [28 il and, 
especially, to having withdrawn the spacelike data ||33ll34ll from the fit procedure. 

One may also note that the 17 accepted decay partial widths |fT9l , which fully determine 
our symmetry breaking (SU(3)/U(3)/SU(2)) parameters, are all well accepted by the fit. At this 
stage, as in our previous study [3], only the p° e+e^ partial width significantly differs from 
its PDG value. 

One should also note that introducing further cross section data sets has to be accompanied 
by the removal of all accepted partial widths |fT9l derived from - or highly influenced by - 
these additional data sets. It is the reason why the number of fit partial widths decreases from 
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17 to 9 as soon as the data on the e+e {n'^/r])'y annihilations are considered in addition to 

e+e^ — > TT+TT^. 



8 Including The e+e {7^^/v)l Cross Section Data 
8.1 Amplitudes and Cross Sections 

Using the Lagrangians given in Appendices |A] and El one can derive the transition ampli- 
tudes 7* P7. The matrix elements are : 



' T(7* - 




= iY [g c^K^{s) 


-(1 


- C3 






< T(7*- 


^07) 


= %Y \g CsMs) 


-(1 


-C3] 






. T(7* - 


Vsl) 


= lY [g csKsis) 


-(1 


-C3] 




e'""'^k^e,{k)paep{p) 



(39) 



where Y = —aemNc/T^f-K, k is the incoming photon momentum (A;^ = s), p the outgoing 
photon momentum {p^ = 0) and Nc = 3. 
We have defined : 



Kp(s] 



E 



DvXs) 



P = 71, r]o, r]s 



(40) 



in terms of the 7 — transition amplitudes Fy^^y (see Eq. (|26l) ) and of the inverse propagators 
Dvii^)- Dpo(s) is given by Eq. (|24|) . while the other inverse propagators are chosen of the 
form s — rriy + imyTv for the narrow u and (p mesons. The functions Hy^^ which carry the 
dependence upon the isospin breaking angles a{s), P{s) and 7(3) are given by Eqs. (1851 ) and 
refer to physical vector fields. We have also defined the constants : 



1 
6 



X {5za + 1) 



6^3 3zA 



3VQ 3za 



(41) 



which are terms deriving from the Caap Lagrangian (see Eq. (l86l)). Defining the three follow- 
ing functions : 



f R^{s) =Y[gc3 K^{s) - (1 - 03)^,] 

cos 6p —sinOp 

sin^P cos6p 
the e'^e^ P7 cross sections are (P = 7r°, i], rj') : 



g csKsis) - (1 - C3)Ls 
g c^Koi-s) - (1 - C3)Lo 



(42) 



a. 



24 



s — rrip 



RUs) 



3af 



s — nip 



\{g CsKp^s) 



C3)Ll 



(43) 



where mp is the mass of the pseudoscalar meson produced in the annihilation process. Kp 
and Lp for P = rj,-)]' are trivially defined from Kq/s and Lq/s using Eqs. (|42)) . 9p, the 
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pseudoscalar mixing angle, is algebraically related with the breaking parameters x and za 
reminded above (see Eq. (E7) in lO). Unfortunately, there is presently no available cross 
section data on the e+e^ rj'^y channel. One can only use the r]''j branching fraction 
|fT9l as constraint ; however, one sees that the cross sections contain also definite constant 
terms (if fits confirm that C3 7^ 1 is significant) and that these constant terms differ for i] and r]'. 
Therefore, the recommended branching fraction [ T9l [351 should be considered with some care, 
until a consistent analysis of the corresponding cross section becomes possible. 

8.2 The Data Set Submitted To Fit 

There are several data sets on the annihilation processes e^e~ 7r°7 and e^e~ — *• 777 
available since 1999, all collected on VEPP-2M accelerator at Novosibirsk. In our analysis, we 
use all the data points up to y/s = 1.05 GeV. 

CMD-2 has recently published a data set on the final states {7i^/r])'j (with n^/i] — > 77) 
from 600 to 1380 MeV 06I with 6 % systematic error. Previously, the same collaboration has 
published data [36] on the 777 final state covering the same energy range and going through 
the mode r] — > 37r° ; these data sets have 6.1% and 4.1% systematic errors, resp. below and 
above 950 MeV. We also use their former data set [[371 on the 777 final state, with r] tt+tt^tt^, 
having a systematic error of 4.8%. 

On the other hand, the SND Collaboration has also recently published [[38l two different 
data sets for the 777 final state with ~ 4.8 % systematic errors, one with 77 — 37r° from 600 
to 1360 MeV, the other with 77 — > tt+tt^tt^ covering an energy range from 755 to 1055 MeV. 
A sample covering the energy range from 600 to 970 MeV for the 7r° 77 decay mode was 
also published [[39l . Other data sets of 14 energy points between 985 and 1039 MeV were 
also published [[40l with both final states (71^/77)7 (and (tt^///) — > 27) and systematic errors of 
2.5%. 

Altogether, these two Collaborations provide 86 measurement points for the e+e" 7r°7 
cross section and 182 for e+e~ 777 for < 1.05 GeV. These data are highly valuable in 
order to build up and thoroughly check our Extended Model in the anomalous sector. 

In this second step, we consider altogether the three annihilation processes e+e^ vr+vr", 
e+e^ 7r°7 and e+e^ 777. In order to constrain more our isospin symmetry breaking 
mechanism, we still include in the data set to be fit a part of the partial widths used in the 
previous Section and already used in [3]. Only 9 pieces of information are now independent of 
the present data : — > n'^'y, K'^ — > if°7, K"^ — > -ft'^7, rj' lu^, (j) r]''^, r]/r]' 77 and 
finally (p tt+tt^ (in modulus and phase), which clearly carry information not statistically 
related with the cross sections considered. 

Instead, as the recommended values [19J for cu/cp — > e'^e~ are information highly influ- 
enced by the set of processes considered (and by the e'^e~ tt+tt^tt*^ data considered later 
on), it is legitimate to let them free. This statement is a fortiori valid for the partial decay 
widths /uj/(j) (7r°/77)7. 

A last remark : the cu and masses and total widths are extracted from the data sets we 
consider. Therefore, it is certainly legitimate to let them vary. At the very end of our procedure, 
when the data on the e+e^ tt+tt^tt" annihilation process will have been considered, we will 
be in position to propose motivated averaged values for this information. Comparing with the 
results derived using the S-factor technics of the PDG [[T9l would become interesting. 
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Concerning our dealing with systematic errors for the newly introduced cross sections, we 
did not find numerically any need to split up correlated and uncorrelated systematic errors, 
which could have allowed increasing the fit parameter freedom. We therefore have simply 
added in quadrature systematic and statistical errors in order to compute the to be minimized. 

Finally, the e"'"e~ — ^ (7r'^/r])'j cross sections submitted to fit have to be corrected for pho- 
ton VP effects. This is done, as in the previous Section, by using the code provided by F. 
Jegerlehner 



8.3 Analysis of the Fit Results 

The fit has been performed and the fit quality information is reported in the third data 
column of Table [U and some fit parameter information is given in the second data column of 
Tabled 

In Tablefll one sees that the fit quality reached for the e+e" 7r°7 (x^/points ~ 61/86) 
as well as for e+e^ 777 (x^/points ~ 129/182) is very good. 

Fig. [U displays together the fit and data for the annihilation process e+e^ 7r°7. The fit 
description is clearly quite good for these three data sets. Moreover, one should note that the 
fit values for the to and peak locations are well centered compared to the data. The situation 
exhibited by Fig. |2] for the annihilation process e+e" 777 is quite comparable. In the to 
peak region, the large fluctuations in the experimental data prevent to be conclusive about the 
detailed lineshape returned by the fit, however, the mass region is nicely accounted for. 

In Table [H one should also note a significant increase of the contribution of the decay 
modes : it practically doubles its value while the number of data is reduced by a half. A closer 
look at the results shows that the contributions provided by 7 out of the 9 modes sums up 
to only 2.5, while x^(0 ^ v'l) = 5.1 and x^iv ~^ 77) = 6.9 , a 2.3a and a 2. 60" effect 
respectively. One may not worry too much about the 2.3a difference from the recommented 
value for Br(0 r/'7) for reasons already sketched ; however, a 2.6ct for Br{r] — > 77) could 
call for some comments. 

When adding more and more spectra to be fit, the weight of isolated independent partial 
decay widths becomes generally less and less constraining. If needed, one may increase the 
weight of the decay mode of concern inside the full x"^- For instance, if instead of adding x'^^-yj 
to the total x^ one adds 4 x xl^-y-f or 8 x xl^jy the distance to the recommended value becomes 
resp. l.lOcr or 0.72a without a significant change to the "decay mode" contribution to the total 

(it increases by 0.3 unit compared to the datum in Table [T])- As this may look artificial, 
a more "natural" way would be to fix za = [Ik/ fn]'^ = 1-495 ± 0.031 at its (experimental) 
central value and use our model equations 



''?77 



TT 



5-2Z ^ + Z . 

cos dp — v2 xsmtJp 



3 3 . , 

. r^-i '''' 

tanb'p = v2 x 

2Z+1 

to connect with the T(P ^ 77) datum [[T9ll : 

r(P - 77) = ^\Gp,,\' , (P = vr°, V, V') ■ (45) 
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One could also derive the 1 cr upper and lower bounds for x and ^^consistent with the fx / f-w 
and r(P — > 77) data and force the fit to stay within these limito Therefore, the problem 
encountered with r(P — * 77) can easily be accomodated without any trouble. 

On the other hand, it is useful to compare the fit parameter values derived in the present case 
with their analogs in the leftmost data column of Table |2l One first notes that the scale factors 
affecting the e+e^ — > tt+tt^ data sets are left unchanged. The second remark concerns the 
physics parameters ; most of them vary well within the quoted la uncertainties. The uncertainty 
for a is improved and its value is still significantly different from 2. The central value for x 
becomes closer to previous fit results using decay width data only [[5l[T4|. 

One should, however, note that C3 becomes inconsistent with 1 by ~ A.ba. This is clearly 
influenced by the anomalous branching fractions but also by the full cross section lineshapes 
of the anomalous annihilation processes we just considered. Forcing C3 = 1 gives a fit quality 
— » ?7'7) = 383.6/469, almost as expectec^lfl. 

It is worth noticing the uj and parameter values returned by the fit. At this step - before 
using 7r+7r"7r° data - we get = 782.45 ± 0.05 MeV and = 8.63 ± 0.08 MeV on the one 
hand, = 1019.25 ± 0.02 MeV and = 4.19 ± 0.05 MeV on the other hand. As will be 
seen later on, these values are modified while including 3-pion data. 

9 Including The e+e tt+tt Data From KLOE 

The KLOE Collaboration, operating at the 0-factory DA<1>NE, has recently published lfT3l 
the spectrum for the e+e" tt+tt^ form factor. Using the Initial State Radiation (ISR) mecha- 
nism, they produced a spectrum covering the region 0.60 — 0.97 GeV with very small statistical 
errors (~ 0.5 %). The systematic errors, also small, are dominant and have been throroughly 
studied M. 

Several source of systematics are reported as (standard deviation) spectra representing frac- 
tions of the measured spectrum for [^^^(s)^. Their [[T3ll Table 1 thus gives the uncertainty due 
to background subtraction (ei(s)), their Table 2 displays the uncertainty due to acceptance cor- 
rections (£2(5)) ; Table 3 and 4 respectively give the error due to detector resolution (£3(5)) 
and the error due to the radiator function effects (£4(5)). From their Table 5, one can derive by 
adding in quadrature the various source of systematics (other than the ones just listed) a global 
scale uncertainty (eo(s)) of 0.76% for \Ft,{s)\'^. 

All these sources of uncertainty should be considered correlated except - maybe - for the 
error due to detector resolution (63) which could have to be treated as uncorrected [42] . This 
structure of systematic errors is clearly complicated and the question is how to deal with the 
error functions ea{s) (a = 0, ■ ■ ■ 4) just defined. 

The most appropriate way seems to follow the method presented in Section [6l i.e. each 
function e^is) is viewed as a gaussian random variable of zero mean and having s-dependent 
standard deviations (named r]a{s)) given by the numbers in the Tables 1-5 of [13J for each 
As bin. More precisely, one may assume < ea{s) >= 0, < [ea(s)]^ >= [rjais]]'^ and < 
ea(si)e/3(sj) >~ SafsSij. This means that these five ea{s) functions play as s-dependent scale 

^'*Unless otherwise stated, the normal running conditions for our fits do not impose bounds to any parameter. 
'^The difference with — 369.5 reported in Table[T]is found at 14.1 while one expects 19.4 ; this shows 
that the minimum is indeed close to parabolic. 
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Parameter 


(NSK Only) 


e^e — ^ tt'^tt 
(NSK Only) 

+ e'^e" — > {Tr'^/rj)^ 


e^e — ^ tt'^tt 
(NSK+ KLOE) 

+ e'^e~ — > {ir'^/r])^ 


All e^e — ^ tt'^tt 
+ e'^e~ — > (7r°/77)7 
+ e'^e~ — > 7r°7r''"7r'~ 


Scale New Timelike 


0.995 ± 0.004 


0.996 ± 0.004 


0.991 ± 0.004 


0.991 ± 0.004 


Scale Old Timelike 


1.007 ±0.009 


1.007 ±0.009 


1.007 ±0.009 


1.007 ±0.009 


Scale KLOE 


— 


— 


1.615 ±0.816 


1.621 ±0.813 


Scale KLOE 1 


— 


— 


-0.041 ± 0.023 


-0.041 ± 0.023 


Scale KLOE 2 


— 


— 


-0.070 ±0.017 


-0.070 ±0.017 


Scale KLOE 3 


— 


— 


0.003 ± 0.006 


0.003 ± 0.006 


Scale KLOE 4 


— 


— 


-0.011 ±0.014 


-0.011 ±0.014 


a 


2.399 ±0.022 


2.356 ± 0.012 


2.361 ±0.011 


2.365 ±0.011 


9 


5.468 ±0.021 


5.574 ±0.019 


5.567 ±0.013 


5.568 ±0.011 


C3 


1.018 ±0.017 


0.943 ±0.013 


0.927 ±0.013 


0.930 ±0.011 


X 


0.935 ± 0.014 


0.904 ±0.014 


0.915 ±0.014 


0.914 ±0.014 


za 


1.577 ±0.020 


1.467 ±0.034 


1.503 ±0.020 


1.496 ±0.018 


zv 


1.509 ±0.020 


1.425 ± 0.045 


1.501 ± 0.030 


1.503 ±0.028 


zt 


1.275 ±0.053 


1.301 ±0.058 


1.340 ± 0.059 


1.332 ±0.058 



Table 2: Parameter values in fits performed including photon VP. The data subsamples in- 
cluded in the full data sample submitted to fit are indicated on top of the Table. The number 
of independent decay widths added to the data sample is 17 (first data column) or 9 (all other 
data columns). For the first 3 rescaling coefficients given in the Table, the corrections are the 
departure from 1, for all others, the rescaling are departures from 0. 

uncertainties. Then a predicted value |F7r(si) p should be associated with a datum modified 
in the following way : 



m,- — > m,- = m,- 



1 ± Yl ^^-{si) 

a=0,l,2,3,4 



(46) 



where each 5Xa{si) is one sampling of the corresponding random variable ea{s) to be fitted. 
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Under these assumptions, the error covariance matrix writes in the usual wajo : 

V^, = + [vM)\FM)?f] + E [Va{s^)\F^{s^)\'] [va{Sj)\F^{s,)\^] (47) 

0=0,1,2,4 

to be inverted numerically for estimation, i, j are bin indices, while is the reported 
statistical error on mj. 

In practice, this turns out to compare the predicted function | (s) p following from a model 
to the modified data : 



1 + Yj QaVa{Si 
0=0,1,2,3,4 



■n (48) 



where the five constants are to be fit. Then, the partial associated with KLOE data set is 

(/i = |F^(s,)|2): 

i,j a=0,--- 4 

In this way, one can check the consistency of with respect to expectations as outlined in 
Section |6] and correct, if needed. 

Finally, as the data under examination have not been corrected from photon VP effects, our 
fitting function is defined by Eqs. (|36l ) and dTT] ) given above. 

We have submitted to fit the data set consisting of all the previously defined data subsets 
plus the KLOE data. The main results are reported in Table [T] (fourth data column) and Tabled 
(third data column). 

One first remarks from Table [H that the fit probability is quite favorable (87 %). The 
contribution from decays is as already reported in Section[8]and calls for the same comment (we 
have x^iv ^ 77) = 2.7 and v'l) = 6.7). Otherwise, the fit quality of the previously 

introduced (Novosibirsk) tt+tt^ timelike data is marginally degraded, while the description 
of the (n'^/r])'j data is unchanged. The contribution from the KLOE data set may look 
large but is still considered acceptable [|42l ; this may reflect the unusual property of being 
highly dominated by systematic errors, always harder to estimate very precisely than statistical 
uncertainties. However, even if the Xkloe large, the global fit does not degrade significantly 
the fit quality of the other data sets and, moreover, the expected physics parameter values are 
not spoiled. 

Some results referring to parameter values are reported in Tabled One observes a change 
in the scale factor of the new timelike data which is shifted by 2a from its expected value 
(0.4 %), while the scale factor associated with the old timelike data is unchanged (~ 0.7a 
from expectation). Out of the 5 scale parameters qi specific to the KLOE data set, we find a 
global scale (go) correction at ~ (1.6 ± 0.8)a from the expected a = 0.76% ; the parameter g2 
associated with the acceptance correction uncertainties is also significantly non-zero while all 
other corrections are small enough to be neglected. 

Parameter values are slightly changed with using KLOE data. One should note the in- 
creased precision on g, the universal vector meson coupling. Finally, Fig. [3] shows superim- 
posed the KLOE data, the fit function, together with the residuals. One clearly sees that the 



This expression corresponds to having treated 63 uncorrelated. In this case, the sum on a in Eq. |49]does not 
extend to a = 3. We tried both possibihties (uncorrelated and correlated) without getting significant differences. 
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description is reasonable. It should be noted however that the worst residuals are in the region 
covered by the £3 correction, i.e. the reported uncertainy on the detector resolution [fT3l . As far 
as we know, this is the first published fit to the available KLOE data. 



10 Including e+e n^n Data 

In contrast with the annihilation channels examined so far, the e+e" vt+tt^tt^ cross sec- 
tion involves the box anomaly sector presented in Section[2l It also involves the non-anomalous 
sector [3J and the triangle anomalies, as the 7r°7 and 777 final states. This process allows to val- 
idate all sectors of the Extended HLS Lagrangian model, as defined at the beginning of this 
paper. 

Without going into much details, one can list the Lagrangian pieces involved : 

• The VVP piece (see Eq. dH)) expressed in terms of ideal vector fields has been already 
given expanded in [5]. This provides diagrams where the photon transforms to physical 
neutral vector fields with amplitudes given in Section [5l The transitions to pn are given 
by: 

SttVtt' (50) 



Cvvp — Ce^^'^^F^yap , C — — ^ ^ C3 



where the physical meson fields uj, p and appear when using Eqs. (l23l) . The 
particles decay with vertices as given in Eq. (1791) . while the neutral physical fields decay 
as follows : 

Cypp = !^pO . ^- g ^+ ^ !^ [^0 _ + . ^- g ^+ (51) 

This provides contributions symmetric in all pvr charge combinations. The argument of 
a{s), P{s) and 7(5) is always the incoming photon 4-momentum squared. 

The same VVP piece provides also a diagram without symmetric partners, where the 
7 — p° transition connects to cjTr^ and the u meson decays in accordance with Eq. (|5T1) . 
In this case the breaking function is «(s+_). The same piece provides also an unusual 
7— p° term connecting with a p° — tt+tt^ in the final state, with the same weight a 



The APPP piece (see Eq. Q) provides a single term 7 — > Stt : 



C^PPP = ieDe>""'f^A^d^7i°d^7i+df37r- , D 



" 4*^^i - C2 + C3J 



(52) 



to the amplitude, and depends on the newly introduced parameter ci — C2. 
Finally, the VPPP piece contributes through the following piece of Cyppp 

Lyppp = lEe^^^'^^Ujld^'K^d^TX^d^'K+dpTl- + ■■■ , E = -TTT^iCl - C2 - C3) (53) 
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limiting oneself to the pion sector for the moment. Using Eqs. (l23l) . one sees that each 
of the 3 physical neutral vector mesons may have a direct 3-pion coupling, however 
modulated by 1, a{s) and 7(5) for resp. the cu, p° and mesons, s is still the incoming 
photon 4-momentum. 



10.1 Matrix Element And Cross Section 

The amplitude for the 7* tt+tt^tt^ transition can be written : 



(54) 



where £^{k) (k'^ = s) is the (heavy) photon polarization vector. Tgym is the symmetric part of 
the amplitude (in terms of the pn 'final' states), while Tp breaks this symrnetry. 
Using the 7 — transition amplitudes given in Section [51 let us define 



N2{s) 



1 



+ ais, 



+ 



1 



+ 



1 



Dpo{s+^) Dp+{s+o) Dp-{s^o) 
1 1 



(55) 



as. 



where the meaning of s+_, s+o and s_o is obvious. Then, we have : 



Us) 



"12^ 
127r2/| 
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1 - -(ci - C2 + C3) + -m'gc^N^{s)N2{s) - -g{ci - C2 - c^)N^{s) 



-m'^gcs 



Dpo{s) 



(56) 

The p propagators have already been defined by Eqs. (|24l) and (|25l) . In order to study 
the three pion final state, we assume Srn^ (see Eq. (|24)) ) to vanish ; the sensitivity to a non- 
vanishing Srn^ is certainly quite marginal as long as r spectra are not considered. 

As in Section [8l we approximate the cu and inverse propagators by : 



DM' 



D4,{q^) =q'^ -ml + im^T^ 
At the chiral point (s+_ = s+o = S-o = s = 0), Tp vanishes and one gets : 

ieNc 



(57) 



-^sj/m (0) 



(58) 



'^If the coupling of the charged p to a pion pair differs from the neutral one, g becoming g + Sg, the last two 
terms in iV2 (s) should be affected by a factor of the form 1 + 6g/g. A test for a non-zero Sg was performed in lO 
using T data and was not found significant. 
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from having required C3 = C4, and having assumed [D^{s = 0))]^^ = [Dp{s = 0))]^^ = — m^. 
Therefore, one recovers the usual WZW term automatically at the chiral limit. 
On the other hand, the differential cross section writes : 



;s'G{x,y)\T,ym + T,\^ (59) 



dx dy 1927r' 

when using the (x and y) parametrization proposed by (43) and summarized in Appendix O 
One should note the abbreviated notation for the functions G{x, y), A^2(s) and N^{s) which are 
actually all functions of x, y and s. 



10.2 Properties of the Matrix Element 

The two pieces Tgym and Tp of our amplitude are given in Eqs. (l56l) . It is easy to check 
that, in the limit of Isospin Symmetry conservation, the non-symmetric term Tp vanishes and 
that the symmetric part named Tsym reduces to only one term (with an intermediate to). 

As vector meson mixing occurs, the symmetric part Tgym clearly exhibits the 3 possible 
7 — y transitions allowed by our Extended Lagrangian Model. The second term of Tgym is the 
more usual one, as it describes the sequence e+e^ followed by 7r+7r~7r° through 

each of the y^pvr possible couplings. One remarks that intermediate = p'^ or: = (p are 
generated by Isospin Symmetry breaking. 

The first term in Tgym is a non-resonant contribution 7 vr+vr"??*^, specific of the HLS 
Lagrangian. It plays, for the 3-pions decay amplitude a role similar to the HLS a (~ 2.4) 
parameter for the e^e^ nn or e+e^ KK amplitudes. The second term is the more usual 
one and should provide the dominant contribution. 

The third term, instead, is the more problematic contact term [43J ; one sees that our model 
predicts that the 3 neutral vector mesons have a direct decay V° tt+tt^tt^ amplitude, with 
the p and contributions weighted by the mixing functions a{s) and 7(5). 

If one follows the educated guess expressed in [2J and thus fixes C3 = C4 = ci — C2 = 1, 
the resonant contact term identically vanishes, while the non-resonant one survives with an 
intensity of -1/2, as supposed in a previous study on box anomalies [|T4l . As we do not expect 
large departures from the quoted guess, one can assert that the resonant contact term should 
indeed be very small. This is confirmed by the data and examined below. 

Finally, one should note the dependence on a(s+_) of the non-symmetric part of the am- 
plitude Tp. Actually, it is this special dependence which defines the non-symmetric amplitude. 
Other approaches generally consider only an oj term, while we also get a p term. Taking into 
account the large width of the p meson, one may guess that this term provides a tiny contribu- 
tion. 



10.3 The Data And The Fit Procedure 

Here also a large number of data sets is available. The most significant are coming from 
CMD-2 and SND Collaborations running on the VEPP-2M machine at Novosibirsk. 

CMD-2 has published several data sets with quoted uncertainties on the cross section merg- 
ing statistical errors and uncorrelated systematic errors. The correlated systematic error, which 
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reflects the uncertainty on the global scale of the cross section, is given separately for each 
CMD-2 data set. The CMD-2 data sets to be considered are : 

• The data set in (23] which covers the cu region and is claimed being affected by a global 
scale uncertainty of 1.3%, 

• The data set given in Il44l which covers the region, with a reported scale error of 4.6%, 

• The (f) region is also explored in [|45l . with a better reported scale uncertainy : 1.9%, 

• The most recent CMD-2 data set is published in ||46l and still covers the region with a 
favorable scale uncertainty ( 2.5%). 

On the other hand, the SND Collaboration has published two spectra covering altogether 
the region from 0.440 to 1.38 GeV. We have considered all data points with y/s < 1.05 GeV. 
These are : 

• Below 980 MeV, the data set in [|471, 

• Above 980 MeV, the data set in [gl 

SND preferred providing the statistical and systematic errors separately in their Tables. 
They claim that their detection efficiency and their integrated luminosity measurement merged 
together are actually a correlated systematic error. Therefore, all systematics, but these con- 
tributions, were added in quadrature to their statistical errors and we treated the rest as an 
uncertainty on the global scale of the cross section. This amounts to 3.4 % below 980 MeV 
[|47]l and we took 5% above [|48]| . 

For both CMD-2 and SND data, we treated the scale uncertainties as explained in Section 
m As in this case, we added penalty terms to the global with the variance just quoted for 
each data set. 

Besides these two groups of data samples, older data sets are worth considering : 

• The former data sample collected by the ND Collaboration, with a reported systematic 
error of 10%, as given by the corresponding Table in the Physics Report by Dolinsky et 
al. El. 

• A small CMD data sample ll50l - 5 measurements - covering the region in between the 
uj and (j) peaks. The reported systematic error is 15 %. 

The former covers the whole energy region between 0.75 and 1.38 GeV, but for our use, it 
was truncated at 1.05 GeV. All data points in this sample lay outside the uj and (p peak regions. 
Therefore, this allows us to constrain the most possible the non-resonant region which is well 
covered by SND, but more poorly by CMD-2. The data points being rather unprecise, there 
was no point in splitting up the systematic error into correlated and uncorrelated pieces ; we 
simply added them in quadrature with the reported statistical errors. We did alike for the CMD 
measurements [i50il . 

Let us mention a DM1 data sample [14911 which covers the energy region between 0.75 and 
1 .098 GeV. It could have been used, however, we are not sure of the systematics all along the 
spectrum. Therefore, we have preferred leaving it aside. 
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We still perform a global fit and, therefore, the data set submitted to fit always includes, 
together with specified 3 pion data subsets, all reported tt+tt^ Novosibirsk and KLOE data 
samples and all e^e^ i'^'^/v)l data sets. 

From the point of view of a fit procedure, fitting 3-pion cross sections using the am- 
plitude given in Eqs. (I56l) poses a problem of numerical analysis. Indeed, as clear from these 
equations, we have to integrate on the Kuraev-Siligadze [43] x and y variables at each step of 
the procedure a function which depends also on the parameters subject to minimization. This 
makes the computing time0 rather prohibitive. 

However, having assumed that charged and neutral p have the same mass and the same 
coupling, the single parameter really influenced by 3-pion data is ci — C2 (see Eq. (l56l) ) ; all 
other parameters in functions to be integrated have already very precise fit values provided by 
the other data sets and are marginally influenced by the 3-pion data. 

Therefore, the method we followed is to tabulate, in bins of A(-yi) =0.1 MeV, the integrals 
over X and y of the G{x, y) function alone and combined with the appropriate products of the 
iVj (see Eqs. (1581) ) and N* functions computed first from a minimization solution without 3- 
pion data. Then one performs the global fit - including 3-pion data - using these functions. 
The fit result is then used to improve the integral calculations and restart a new minimization 
step. A priori, this procedure has to be repeated until some converge criterium is fulfilled. 

We expected convergence in a very few steps. However, it so happens that there was no need 
to iterate the procedure, as the was never changed by more than ~ 0.5 unit, compared to 
some 'exact' computing method, where the integrals were computed at each minimization step. 
As the fit parameter values were not found to get significant corrections, we found justified to 
use this simplified method for performing our fits. 

10.4 Exploratory Fits 

Up to now, we have always used all available data sets for the various annihilation channels 
we have examined. The 3-pion data faces us with having to perform a choice between data sets 
which has to be motivated. It is the reason why some exploratory analysis of the 3-pion data 
sets has been worth performing. 

We start by considering aU CMD-2 [I23lll3ll5l|46l and SND [HTlllSl data sets in the less 
constraining pattern of our fit procedure ; namely, together with only the so-called old and new 
Timelike data sets collecting the standard e+e^ tt+tt^ data sets. In this case, we have to use 
ll3l also the full set of 17 accepted decay modes lfT9l . Indeed, doing this way : 

• We lessen at most effects of e+e^ (n^ /r])'j on the value of C3, which plays some role 
in the 3-pion cross section, 

• We avoid being sensitive to the KLOE data and their particular systematic errors. 

'^In order to analyze the behavior of each parameter or subset configuration, we need to run the code a large 
number of times. Without 3-pion fitting, such a run lasts about 2 minutes. Including the 3-pion samples, the 
fastest 'exact' method we found provides a waiting time of about an hour The simplified method sketched below 
reduces the computing time to about 5 minutes per run. 
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corr (C3, Cl — c2 ) 


0.027 


0.252 


0.523 


0.152 


0.473 


xVdof 
Probability 


454.73/413 
7.7 % 


515.48/564 
92.9 % 


815.40/735 
2.1 % 


667.21/646 
27.4 % 


637.90/653 
65.6 % 



Table 3: Fit results using three-pion data sets within global fits. In the first line the number 
of decay widths is indicated in boldface. The so-called 'conditions' are the correlated scale 
uncertainties which are considered as data and fit. Comments are given inside the text. The 
partial in lines flagged by * do not include the additional contribution of the scale factor. 
These are, nevertheless, counted inside the final together with the number of conditions. 



While comparing with other models and fits, one should keep in mind that, even if 
highly relaxed, our model is, nevertheless, sharply constrained. It is actually its lack of flexi- 
bility which imposes some good enough understanding of the experimental uncertainties and. 
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on the other hand, discloses the correlations between the various data sets imposed by the 
(common) underlying physics. 

The fit pattern just described leads to the results reported in the first data column of Table 
|3l The fit qualities for the CMD-2 data sets given in [|44ll46ll are good (x^/npoints ~ 1) while 
the fit qualities of the data sets given in ||23l |45l are worse (x^/npoints ~ 1.5 -h 2). The fit 
quality of the old ND [51] and CMD [|50l data samples is also quite good (x^/npoints ~ 1 
altogether). The low energy SND data set [|47ll returns a good fit quality (x^/npoints ~ 1), 
while the higher energy SND data set ll48l is found worse (x^/npoints ~ 1.5). However, the 
global fit quality (x^/npoints ~ 455/413) corresponds to a 7.7% probability. 

This relatively poor probability leads us to carefully examine the behavior of the CMD-2 
and SND 3-pion data sets separately. 

The next step relies on the following remark. All data sets - except for 3-pion data - 
allow to fix all parameters of our model, but ci — C2. As can be seen from the third data column in 
Table[2l these are already known with an excellent accuracy, especially the parameters expected 
to be influenced by 3-pion data {g, c-s for instance). Looking at Eqs. (l56l) . one can easily see that 
the missing parameter ci — C2 occurs in such a way that it is highly influenced by the invariant 
mass region outside both the to and peaks. Therefore, one can guess that performing the 
global fit with the ND [|5T| and CMD [|50l data samples only should allow a measurement of 
Ci — C2 with already a challenging accuracy. On the other hand, the Breit-Wigner mass and 
width for the to and meson are accurately fit from the {Tc^/r])'y data. 

Therefore, one can construct an accurate 3-pion cross section, valid for the whole energy 
range from threshold to above the (p mass, which feeds in the whole correlated physics infor- 
mation. This cross section can be used as a prediction to be compared to the data sets provided 
by the CMD-2 and SND Collaborations separately. 

The global fit following this pattern has provided results shown in the second data column 
of Table [31 The global fit quality is well illustrated by the fit probability which reaches 93%. 
A closer look at the fit properties of the previous data sets shows that they are still optimally 
accounted for (compare this data column with the third data column in Table [2l). The fit quality 
reached for the ND + CMD data set is also very good : x^/npoints ~ 26/37. 

The accuracy reached for g and cs is unchanged, however this fit provides a value for 
ci — C2 — 1.2 with an accuracy of ^ 4%. Therefore, this allows us to define precisely the 3- 
pion cross section. This function is displayed in Figured together with the CMD-2 and SND 
data superimposed in the to region and in the region. 

The left plots in Figure |4] show the case for CMD-2 data ; one can see that the data are 
already in good agreement with what can be inferred from ND+CMD data together with the 
physics of the other annhilation processes. The peak value and lineshape for the cu region are 
in accord with the prediction ; one may note, however, minor departures far on the wings. 
Concerning the (p region, the lineshape is in reasonable agreeement with expectations, and one 
may infer that a fit will result in a good description. 

Rightmost plots show the case for SND data ; at both the to and peaks, one observes that 
both tops are ~ 10% too large compared with the same expectations as before. However, a 
qualitative observation of the behavior in the region, allows one to infer that one should run 
into difficulties to account for both top and wings with only global (constant) rescalings. 
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This conclusion is substantiated with performing the global fit with all SND and CMD-2 
data, as reported in the third data column in Table [3l In this case, the global fit probability 
reached is ~ 2% ; the scale factors returned, 0.98 ± 0.03 for the low energy data sample 
[|47l and 0.91 ± 0.04 for the higher energy one P8ll are in good agreement with expectations. 
However, this poor fit probability, associated with poorer SND x^'s, indicates that the problem 
is not solved by means of constant rescaling factors. 

Of course, this follows from having requested consistency of the 3-pion data sets, not only 
with each other, but also with all physically related annihilation processes. Indeed, one clearly 
understands that a common fit, using only the CMD-2 and SND 3-pion data, should certainly 
succeed with a solution intermediate between CMD-2 and SND data ; but this success would 
hide the consistency issue with the rest of related processes. 

Therefore, the two global fits performed with merging SND and CMD-2 data tell us that 
there there is some inconsistency between them. Now, it remains to check separately SND 
+ND +CMD data, on the one hand and CMD2 +ND +CMD data, on the other hand, in order to 
make a motivated choice. 

We have performed the global fit with all SND+ND+CMD data ; the results are reported in 
the fourth data column in Tabled The global probability becomes more reasonable (27%) - 
from having removed CMD-2 data - and the fit rescaling factors are 0.96 ± 0.03 and 0.89 ± 
0.04 for resp. the low and high energy regions. The fit quality exhibited for the low energy 
sample may look reasonable, however, the corresponding information for the higher energy 
data sample goes on looking poorer. This should reflect that the 10% scale correction, valid on 
the peak, is not appropriate on the wings. 

The corresponding information from the global fit performed with all CMD-2+ND+CMD 
data is reported in the last data column of the same Table ; one gets a global fit probability of 
66 %. This result, together with what is displayed in Figure IH shows that ND+CMD data, on 
the one hand, and CMD-2 data, on the other hand, are quite consistent with each other. 

Clearly, as the physics contents of the CMD-2 and SND data sets are the same, the issue just 
sketched cannot be attributed to our model. Because of Figure |4] and because of the respective 
probabilities, one has to choose among CMD-2 data and SND data, and the best motivated 
choice is the CMD-2 +ND +CMD solution. This also turns out to remark that Table [3] exhibits 
a poor consistency of SND data compared with ND +CMD. SND data are certainly useful, 
however not in the framework of a global fit which requires a good control of systematics in 
order to get reasonable consistency with all other physics process measurements. 

10.5 The Solution With CMD-2, ND and CMD Data 

Presently, the choice to remove the SND samples for 3-pions data considered in the global 
fit, seems the best motivated one. The main fit results are shown in the rightmost data column 
in Table [3l The fit is shown with data superimposed in Fig. [51 Top plots show a zoom on 
the peak regions, while the downmost plot focuses on the region outside the peaks. The fit is 
clearly satisfactory. 

Let us summarize the other physics results of interest : 

• The rescaling factors for all Novosibirsk and Frascati data reported in Table [2] are recov- 
ered with changes affecting the last digit, i.e. far inside the reported errors. 
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• the value for the HLS parameter a = 2.365 ± 0.105 is also unchanged. 

• The breaking parameters za, zy and zt were also left unchanged by the fit procedure, 
with unchanged uncertainties, compared with the third data column in Table [2l 

• The scale factors found are 0.996 ± 0.012 for the data set in [23J (expected departure 
from 1 : 1.3%), 0.975 ± 0.020 for the data set in B6l quite in accord with the expected 
departure from 1 (2.5%). The other rescaling factors are resp. 0.953 ± 0.035 (data from 
Bi] with expected departure from 1 of 4.6%) and 0.974 ± 0.016 (data from [45J with 
expected departure from 1 of 1.9%). Stated otherwise, the first 3 fit scales are 1 a at most 
from expectations, while the last scale is found at 1.25 a. This also indicates that there is 
no obvious signal of missing variance for the 3-pion CMD-2 data samples. 

The fit quality for the ND data set in [[5T]| and the old CMD data set IfSOl is quite satisfactory. 
The three CMD-2 data sets collected around the peak benefit from a x^/npoints ~ 1 ; 
instead, the for the data set collected at the cu peak has a much poorer quality (x^/npoints = 
25/13). A closer look at the results shows that more than half of the comes from 3 points 
on the wing^. However, as the global lineshape looks well reproduced at both peaks and in 
between them and also because of the good global fit quality, we consider the fit to the CMD-2, 
ND and CMD data as satisfactory. 

Other fit information, reported in the last data column in Tables [U and [2l allows estimating 
the effects of having included 3-pion data inside our fit data sample. Now, with the full in- 
formation displayed. Table 1 clearly illustrates that the only noticeable changes (x^ increased 
by ~ 4 units) concern the 7r°7 and 777 data samples. Table |2] instead shows that there is no 
modification beyond the 1 a level. 

11 A Few Numerical Results 

With the fit presented just above, one ends up introducing data affecting the cu and mesons. 
Our global fit already provides interesting physics results which should not be affected by the 
forthcoming steps of our study. 

11.1 id and (J) Masses And Widths 

We have defined the cu and propagators as fixed width Breit-Wigner expressions (see Eqs. 
(l57l)). Our fits widely illustrate that there is no need to go beyond this approximation. As we 
account for both statistical and systematic errors, the uncertainties we quote fold in both kinds 
of errors and are thus directly comparable to RPP |fT9l information. Moreover, it should be 
noted that the RPP data on u and masses and widths mostly relies on the data we have used. 

The measurements located at ~ 0.76, ~ 0.77 and 0.80 GeV contribute resp. for 5, 4 and 6 units to the . 
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Our final results are : 





= 782.42 ± 0.05 MeV 


, RPP : 




= 782.65 ±0.12 MeV 




= 8.700 ± 0.084 MeV 


, RPP : 




= 8.49 ± 0.08 MeV 




= 1019.173 ±0.015 MeV 


, RPP : 




= 1019.455 ± 0.020 MeV 




= 4.259 ±0.036 MeV 


, RPP : 




= 4.26 ± 0.04 MeV 



(60) 



We thus find the lu mass 2 a below RPP average value and we get an uncertainty twice smaller. 

is found at 2.5cr from its recommended value [T9], but in close agreement with the fit results 
of both the CMD-2 and SND Collaborations reported in the Review of Particle Properties [19J. 

The (j) mass is instead found at 0.282 MeV from its recommended value |fT9l - about 20 a ! 
Interestingly, our extracted (p mass is also significantly smaller compared to the values extracted 
by CMD-2 and SND from the same data [fT9l . This might partly reflect the contributions of 
anomalous terms and the effect of vector meson mixing. Moreover, as parametrization of the 
propagator, we used a fixed width Breit-Wigner function, while CMD-2, for instance, used 
a varying width Breit-Wigner expression. Finally, P^ is found in perfect agreement with the 
RPP for both the central value and its uncertainty. 

One should note that our global fit is expected to : 

• lessen the statistical error. Indeed, each parameter is constrained simultaneously by all 
data sets where it plays a role. This is, of course, true for the mass and width of the uj and 
(j) mesons which sharply influence the 7r°7, 777 and the tt^tt+tt" final state descriptions ; 
moreover, the sharp drop observed in e+e^ tt+tt^ is also influencing the a; meson 
parameters. 

• perform the folding of systematic and statistical errors in accord with the whole knowl- 
edge of these uncertainties provided by the various experiments. 

Therefore, the net expected result is an optimal folding in of all reported sources of errors, even 
when data are coming from different triggers and/or groups. 

11.2 Contact Terms And Anomalous Lagrangians 

A priori, the Extended HLS Model we use has two kinds of contact terms. One of the form 
'jPPP describes the residual coupling of photons to pseudoscalar meson triplets, the other of 
the form VPPP the direct coupling of a vector meson to a pseudoscalar meson triplet. The 
former is given (see Eq. (|56l) ) by : 

= 1- ^(ci -C2 ±03) = -0.61 ±0.10 , (61) 

and the latter by the coefficient : 

= -^(ci -C2 -ca) = -0.63 ±0.31 ; (62) 
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this second coefficient is additionally modulated by 1, a{s), 7(5) for resp. the u, p, (j) mesons 
when coupling to a pion triplet. The numerical values just given are derived from the infor- 
mation in Tableland include the effects of the covariance term < 5{ci — C2) 80^ >. These 
numbers are very close to the educated guess in ^ {C^ = — 1/2 and Cy = 0). 

One can conclude that the contact term is certainly significant, while Cy has only a 2cr 
significance. 

In view of these results, one is interested in seeing what happens to the global fit, if one 
fixes C3 = Ci — C2 = 1 in the model. The result is shown in the second data column of Table IH 
The subsamples not shown keep, roughly speaking, their usual contributions. In this case, as 
could be expected, the fits to the e+e^ 7r°7 and e+e^ — > 777 cross sections are significantly 
degraded and the global fit probability allows to reject this solution. 

If one fixes Ci — C2 = 1 and leaves free C3 (third data column in Table HJ, the good description 
of 7r°7 and 777 data sets is recovered as could be expected, but this is done at the expense of a 
worse decription of tt^tt+tt^ data. The situation is similar if, instead, one fixes C3 = 1 and let 
free ci — C2. 

One may conclude from the results collected in Tabled that it is highly meaningful to let 
ci — C2 and C3 vary. This means that the violation of Vector Meson Dominance affects all sectors 
of the anomalous HLS Lagrangian [|71 |2l ; this violation is only weaker in the Triangle sector 
than in the Box sector. 



Global fit 
(X^ contributions) 


C3, ci - C2 free 


C3, ci — C2 fixed 


ci — C2 fixed 


C3 fixed 


7r07 (86) 


65.66 


106.80 


80.72 


69.41 


777(182) 


135.20 


156.14 


136.57 


162.15 


ttOtt+tt- (130) 


137.06 


205.32 


171.24 


159.52 


Probability 


65.66 % 


0.01 % 


12.5% 


12.7% 



Table 4: Influence of the parameters C3 and ci — C2 on the fit results. Each entry displays the 
contribution to the total of the quoted subsample. Last line provides the global fit probability. 
For the first column C3 and ci — C2 values are given in the last column of Table [3] ; the second 
column gives the best solution for C3 = Ci — C2 = 1. In the third data column one has set 
ci — C2 = 1 and the fit has returned C3 = 0.830 ± 0.002. In the fourth data column one has fixed 
C3 = 1 and fit ci - C2 = 1.382 ± 0.038. 
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11.3 The Mixing 'Angles' 



Most of the present data sets examined altogether are dominated by the u and resonances. 
Using cross sections instead of only partial decay widths may produce changes in the mixing 
functions a{s), P{s) and 7(5) which, as noted before, can be interpreted as resp. the co — p, 
(j) — p and — mixing 'angles'. These are shown in resp. Figures [6l|7] and [H 

Compared to using only partial width decays together with e+e" tt+tt^ cross section 
data, one observes that the functions a{s) and are unchanged (compare to Figure 7 in [3]). 
The function 7(5), instead, gets significant changes. Figure[8]shows that the — mixing angle 
becomes positive for s > 0, while it was negative when using partial width decays. However, 
the variation of 7(5) remains significant between the uj and cp regions. The values which can be 
read off Figure [8] give 7 ~ 2.75° at the uj mass and 7 ~ 3.84° at the cp mass. As, now, one relies 
on the largest possible data set, this result should be considered as superseeding the function 
l{s) in 

12 Validation Of The Model : The r//r/' tt+tt 7 Spectra 

The amplitudes for 77/77' 7r+7r^7 at the chiral point can be derived from the WZW 
Lagrangian [fTOl [TTll . However, as there is no theoretical knowledge of their momentum de- 
pendence, modeling the behavior outside the chiral point is unavoidable. As shown by Eqs. 
(fT9l ). the Extended HLS Model allows to recover the chiral limit and the question is whether 
the predicted spectra and partial widths are in agreement with data. Such a study was already 
performed a few years ago [14J assuming, as inferred by |l7l|2l, that ci — C2 = C3 = C4 = 1. 

At the point where we are, all parameters of our Extended HLS Model have definite values 
and, thus, all vector meson and contact term contributions have no longer any free parameter. 
The question is now whether the spectra and the partial widths, which can be algebraically 
derived, fit the known physics information. 

The most reliable physics information for these 2 decay modes is certainly their partial 
widths [fT9ll . There exists two spectra giving the photon momentum distribution in the rj rest 
frame ||52l[53ll ; however, the relevant information should be read off plots giving distributions 
for the measured spectrum and for the acceptance/efficiency function. 

Understanding the dipion invariant mass spectrum in the rj' 7r+7r^7 has been addressed 
several times, as most groups were claiming that the p mass peak was observed shifted. Most of 
these spectra were published only as figures [ f54l[55l[56l [571 [581 ; only Crystal Barrel provided 
the data points, feeding in the information concerning acceptance and efficiency corrections 
|[59l ; some other spectra were only published as preprints or PhD theses. All these spectra 
were discussed in [il4il and shown to exhibit very different qualities due to statistics fi54il or, 
sometimes, to obvious biases |[57l . etc... 

The issue was whether the p peak in the rj' 7r+7r^7 spectrum is shifted. Crystal Barrel 
|[59l clearly proved that the p peak location was indeed shifted by ~ 20 MeV. The study in 
|[T4l later proved that this shift was actually the way found by fit procedures to account for 
a missing constant term which results in a distorsion of the p lineshape. This distorsion is 
produced by the contact term 'jir'^Tc'r]' which adds up with p meson contribution. Accounting 
for the contact term, [[T4|| proved that the p peak in the 77' spectrum is at the location expected 
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from e+e tt+tt data. It is this question which is revisited within a framework where the 
condition ci — C2 = C3 = C4 = 1 is violated and the vector meson mixing at work. 



12.1 Amplitudes For The 77/77' ^ tt+tt 7 Decays 

Using the Lagrangian pieces given in Section [3l it is now easy to compute the transition 
amplitudes involved in the rf/rj' ^ 7r+7r^7 decays. These can be derived from : 

T(r/8 ^ 7r+7r-7) = A Fsis)e'^''''''e^i^)q,p-p+ 



T(r/o ^ 7r+7r-7) = xV2 A Fo(s)e'^'^-^e^(7)g,p-p+ 
Having assumed C3 = C4, the octet and singlet functions write : 



(63) 



A = 



1 



2 

ieNr 



(64) 



where q is the photon 4-momentum and s = (p+ ) ^ is the dipion invariant mass. Neglecting 
terms of order greater than 1 in the mixing parameters a, (3 and 7, one has : 



as 



1 + 



as 
1 + ^ + 



2^/2 
V2 



a(s) 1 ^ 2^2 (3{s 



3 D^{s 
a(s) 1 



3za D^{s) 
V2 P{s) 



(65) 



which fulfill : 



cUs) 



za 



3 D^{s) ?>zaD^{s) 
1 1 



DJs) 



(66) 



The (f> contributions in these expressions can be dropped out, because of the phase space 



cuts at s 



m^, or s 



m 



^, much below the mass. Defining : 







f cosOp —sin6p 






\^ sin6'p cos 9p 




the decay amplitudes for 77 and r]' can be written : 

and the decay partial width can be obtained by integrating : 

dT{X 7r+7r"7) 1 aem 



(67) 



(68) 



\Fx{s)Wlpl 



X = 7], 7] 



from the two-pion threshold to the i] or r]' mass. One has also defined q^ 
andp^ = V ~ 4m2/2. 



X 



(69) 

s)/2mx 



34 



12.2 Properties of the Amplitudes 



The amplitudes just given have several interesting properties. As in the model developped 
in [T4J, one finds in Eqs (|64|) a vector meson amplitude which adds up with a non-resonant 
amplitude. 

One might be surprised that these expressions do not depend on ci — C2 as the decay ampli- 
tude for e+e" — > tt^tt+tt^. This is due to an unexpected conspiracy between the VPPP 
and APPP couplings (see Eqs. (fT6l) and (flTl) ) which both play a role in the transitions 

rj/rj' 7r+7r~7. 

As for the amplitude for e+e^ tt+tt^, the dipion invariant mass distribution exhibits 
an interference pattern between the p and uj meson contributions. However, all experimental 
data sets listed above indicate that this usual interference pattern should not be met, as there is 
no observed drop of the distributions around ~ 782 MeV. The single difference between the 
e^e~ — > TT+TT^ and rj' 7r+7r^7 processes is that the term of order a is weighted by 1/3 in the 
•q' decay amplitude. It is thus interesting to see whether this different weighting alone allows to 
cancel out the (expected) drop in the dipion spectrum. 

12.3 Comparison With Data 

In Figured one displays the dipion mass spectra in the rj' 7r+7r~7 decay with, superim- 
posed, the predictions coming out of the model fitting e+e~ tt+tt^, 77^7, 177, tt+tt^tt^. We 
thus show the case for ARGUS data [56J, for an experiment run at Serpukhov on the Lepton F 
facility [|58l with relatively large statistics, and for Crystal Barrel [|59l , which is certainly the 
most precise rj' spectrum. The last plot in Figured is the spectrum of the photon momentum in 
the rj rest frame for the rj 7r~^7r~7 decay [l53l, together with our prediction. 

The agreement is satisfactory. One should note that some effect of the uj — p interference 
on the predicted lineshape can be detected but is not contradicted by the data. Actually, the 
data are binned and the curves shown are not averaged over the bin widths. One clearly sees 
that the peak location, as well as the global lineshapes are well predicted by our model. One 
may remark that the 20 MeV shift of the p peak confirmed by [59] is correctly reproduced. 
Therefore, having C3 7^ 1 and vector meson mixing does not degrade the agreement already 
obtained in lfT4l with a much simpler model. 

However, the quality of these data is not good enough to allow including the corresponding 
spectra into the global fi0 

Other pieces of information are the partial widths for 77 and rj' decays into this final state. 
The predicted values for these are : 



Actually, the Crystal Barrel spectrum could be safely fit, except for a single point - at 8 12.5 MeV - which de- 
grades severely the without changing the parameter values at minimum. The Crystal Barrel spectrum is actually 
the merging of 4 spectra collected in pp annihilations at rest. Two of these (frompp uj{^ tt'^'j /tt'^tt^tt^ )'/]') 
exhibit this faulty point, for two others {pp ■n^ii^rj' and pp Tr+Tr^ry') the corresponding measurement is 
located on the predicted curve. We have preferred avoiding to include a truncated spectrum inside our data set. 




TT 



7r+7r-7) = 53.11 ± 1.47 keV 



+7r-7) = 55.82 ± 0.83 eV 



(70) 
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quite close to the corresponding recommended values fl9\ : 60 ± 5 keV and 60 ± 4 eV. The 
predictions are clearly quite accurate and provide support to the parameter values we derived 
from our global fit. This also indicates that we have exhausted our parameter freedom with 
the data set we already considered. Indeed, including these two pieces of information inside 
our global fit does not return results substantially different from the predicted values given just 
above. Therefore, the rj/rj' ^ 7r+7r^7 decay modes represent a good test and better data on 
these decays would certainly be valuable in order to further constrain our model. 



13 Decay Partial Widths Of Vector Mesons 

Having allowed C3, ci — C2 7^ 1 in the HLS model [|2l, one allows other anomalous con- 
tributions than VVP to partial widths. Then, beside the diagrams represented by Fig. (flOk ). 
radiative decays of vector mesons get also contributions from diagrams as shown in Fig. (flOb). 
These contributions are proportional to 1 — C3 and are generated by the AAP Lagrangian piece 
of our Extended Model. Correspondingly, the decay widths to 3 pions get their dominant con- 
tributions from diagrams generically represented by Fig. (fTOb). but also subleading V PPP 
contributions proportional to ci — C2 — C3 from diagrams shown in Fig. (flOb ) and APPP 
contributions (see Fig. (flOb)) proportional to 1 — 3(ci — C2 + C3)/4. 



13.1 Radiative Decays Of Vector Mesons 



As far as radiative decays are concerned, let us denote gvP'y the couplings generated by the 
VVP couplings and given in Appendix E of [3J where C3 = 1 was assumed. In our Extended 
HLS Model, they become : 

Gvp, = 9VP, C3 , V = K*', K*^,p^ . (71) 

However, the couplings to ideal fields (Vj = pj, uj, (f)j) yield additional corrections : 



eDFv^{ml) 

dVinO^ C3 : ^ , V 



Gvjr^s-r — gViri«'^ C3 



Gvirf-y — gVirjO^y C3 — 



3 
eP 

e X D 
"9 



mr 



za 



(72) 
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5ZA+ 1 
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The dominant VVP contributions are first order in g and are weighted by C3 ~ 1. The 
additional AAP contributions are of order and are additionally suppressed by (1 — C3). 
In the annihilation processes e+e^ 7P, the AAP terms are subleading and are actually 
absorbed by the intermediate photon vacuum polarization in the e+e^ annihilation amplitude. 
From a numerical point of view, these additional terms, all proportional to the 7 — transition 
amplitudes Fy^(my), give a negligible contribution. These transitions amplitudes are given in 
Section [51 In order to compute the radiative decays of vector mesons, ideal field combinations 
have to be constructed [3J, as reminded in Section HI Also, pseudoscalar singlet and octet 
couplings have to be combined in order to derive the r] and rj' couplings [3]. The dependence 
upon the mixing 'angles' of vector mesons is hidden inside the gvp-y and Fy^ functions. 
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13.2 Three Pion Decays of Vector Mesons 



The coupling constants for three pion decays of neutral vector mesons can be derived in 
close correspondence with the 7* tt+tt^tt^ amplitude constructed in Section [lOl Let us 
define : 

3 9 

M{x, y, s) = -m'^gcsNiix, y, s) - -g{ci - C2 - C3) 



C 



E 



1 - ^(Cl -C2 + C3) 



(73) 



fully displaying the dependence upon the Kuraev-Siligadze f43l variables. The function N2{x, y, s) 
has been defined in Section [10] and = ag'^f^. Then, the amplitudes for V 7r+7r~7r° are, 
at leading order in the symmetry breaking parameters : 



Auj{x,y) 
A<p{x,y) 



E 



E 



M{x,y,ml) + e^C 



mf: 



Ap{x,y) = E 



7 m , M X, y, + e C 

^ ^ mi 



a{ml)M{x, y, ml) + -m^gc2,N'i{x, y, ml) + e^C 



(74) 



Fp~m) 

mi 



where N-i{x,y, s) has also been defined in Section [TOl and hides a dependence upon a(s+_) 
(s+_ being the tt+tt^ squared invariant mass of the decay products). One may observe that 
the dominant terms for p and (p are generated by vector meson mixing and that there is a 
specific term in the p amplitude compared to those for uj and </>. One should also note that the 
APPP coupling generates a subleading term (of order e^) which exists even if there were no 
vector meson mixing. The present model is more complicated than the one developped by ||43l 
because of these APPP terms and of the vector meson mixing. 

Finally, in terms of the amplitudes given in Eqs. (l74l) . the 3-pion partial widths are given 
byHH : 



mi 



Trim 



7687r3 



dx dy G{x,y)\Av{x,y)['' 



(75) 



where the function G{x, y) [|43l has been reminded in Appendix O 



13.3 Effects Of A 0{p'^) Lagrangian Piece 

In [[6OI . Harada and Yamawaki have studied the Wilsonian matching of the HLS model with 
QCD. For this purpose they have identified 35 0{p'^) Lagrangian pieces of the HLS Lagrangian 
which are provided explicitly in [2J. Among these, one is of special interest for our purpose, 
the so-called 23 term which is of concern for the 7 — V and W — V transitions. Discarding 
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effects of SU(3) symmetry breaking, the additional Lagrangian piece writes : 

< 

>Cf ^ = Z3gg2Vud d^W^{df,p^ - d^p;,) + ■■■ 
= zzgg2Vud d^W-{d^pt - d^p^) + ■■■ 

(76) 

using notations already defined in |l3l and reminded above. This turns out to modify in exactly 
the same way the transition amplitudes for 'j — V and — V by adding a s-dependent term to 
the constant part of the couplings. More precisely, one gets (see Eqs. (30) in [3J for r decays 
and Section [5] above for e+e" annihilations) : 

where Z3 can be related with the Low Energy Constant (LEG) Lg of the Chiral Perturbation 
Theory. When on mass shell, the additional factor is simply (1 — g'^z^) and can be guessed 
p] of the order 1.10. Indeed, even if z^ is expected small, the g'^ term produces an important 
enhancement factor. 

However, in processes like e~^e~ annihilations or r decay, /J or fv-y always come in com- 
bination with the loop dressing functions [[31 Hwis) or ny^(s) which contain subtraction poly- 
nomials to be fitted. Therefore, the 23 contribution comes in entangled with the first degree 
term of the subtraction polynomial and cannot be singled out. As a consequence, the fit value 
for g may absorb effects of the z^ correction and may produce too low values by as much as 
10% for the amplitudes [H. 
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In the e~^e^ annihilation channels considered in the present work, our fit procedure is 
sensitive to the product Fy-y x g for each resonance rather than to the Fy^'s and g separately Ell- 
The fit quality we have reached allows us to consider that these products are well understood. 
Therefore, all pieces of information relying on these products can be considered reliable. This 
covers all products of widths like T{V P7)T{V e+e^) or T{V 7r+7r-7r°)r(y 
e'^e^). Ratios of widths performed from the T{V P7)'s and T{V 7r+7r^7r°)'s are also 
free from the disease mentioned above and can be considered as secure. It is the reason why 
we will not go beyond these products and ratios until a satisfactory solution to the issue just 
mentioned is found. 

Actually, some possibilities exist to solve this ambiguity, at least in principle. Indeed, 
processes like e+e^ /^^A* are sensitive to the [Fy^J^'s and could be used. Some scarce 



our full data set, only the radiative decays of the K*^, K*^, do not belong to this category, but are not 
sufficient to modify the picture. 
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(accurate) data (3 measurement points) at the mass are reported [62J but, unfortunately, noth- 
ing with the required accuracy is reported around the uj mass. Another possibility could be 
to examine channels like e^e~ uj{-^ 7r°7)7r° which are rather sensitive to some of the 
X (yi^'s. Finally, one could also fix the leptonic widths as coming from the fits of the Par- 
ticle Data Group |fT9l , which basically corresponds to our results published in [|3]|. This last 
solution, even if legitimate, is not completely satisfactory. 

It should be stressed that the effect mentioned is actually common to any parametrization of 
e^e~ and r data and is by no means specific to the HLS Model (see, for instance, the discussion 
in ll43l ). On the other hand, reporting on products and ratios is the standard way experimental 
groups proceed with the kind of data we are dealing with and we may compare to these. 



Data 


Our Fit 


PDG 2008 


r'(p° tttt) xlO^ 


4.72 ±0.02 


4.876 ± 0.023 ± 0.064 


V'{uj tttt) xlO^ 


1.146 ±0.057 


1.225 ±0.058 ±0.041 


r'(pO ^7r07) xlO^ 


1.875 ±0.026 


2.8 ±0.4 


V\lo 7r°7) xlO^ 


6.80 ±0.13 


6.39 ±0.15 ** 


r'(0 ^ 7r°7) xlO^ 


4.29 ±0.11 


3.75 ±0.18 


r'(pO ^ 777) xlO^ 


1.05 ±0.02 


1.42 ±0.10 


T\uj r]-i) xlO* 


4.50 ±0.10 


3.31 ±0.28 ** 


r'(0 r/7) xlO^ 


4.19 ±0.06 


3.87 ±0.07 ** 


r(p° ^ tttttt) xlO^o 


9.03 ±0.76 


45.8t?|^ ± 15.6 ** 


T'{uj tttttt) xlO^ 


6.20 ±0.13 


6.39 ±0.10 


r'(0 tttttt) XlO^ 


4.38 ±0.12 


4.53 ±0.10 



Table 5: Results for T'iV ^ f) = TiV e+e-)T{y f)/Tl^t for each each V vector 
meson and each final state /. The most significant results are flagged by **. 

13.4 Results For Partial Widths 

In Tabled we give our fit results for : 

T'iV ^ /) = viv ^ e+e')v{y ^ /)/rL 

and provide the measurements collected in the RPP |fT9l . The results on which we report have 
been derived from a fit to all e+e annihilation data, except for KLOE data [[T2|. by fixing the 
scale corrections to zero in accordance with Section [6l 
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One thus notes that, for most of the reported products, one gets resuks in accordance with 
what is extracted from the same data. One should remember that the errors we quote fold in 
systematic and statistical errors and can be compared to the averages performed by the Particle 
Data Group. 

Concerning p° Tiim, our result, mostly determined by CMD-2 data, is in agreement with 
the SND result, but with a much lower central value and a much better accuracy. Significant 
differences mostly concerning uj /cp 777 can be read off Table [51 They correspond to per- 
forming a simultaneous fit of all the involved data which happens to be of very high quality as 
shown in Section [10] above. Finally, as a general statement, our uncertainties reported in Table 
[5] are either of magnitude comparable to the existing information or of much better accuracy. 



Data 


Our Fit 


PDG 2008 


T{p^ — >■ '7r7r)/r(p*' vrvrvr) 


(5.26 ±0.45) 10^ 


> 100 


r(a; — s> 7r7r)/r(a; — ^ tttttt) 


(18.5 ±1.0) 10-=^ 


(17.2 ± 1.4) 10-3 


r(|o'^ — >■ 7T'y)/r{p^ — s> tttttt) 


20.91 ± 1.60 


[5.94 ±3.36] ** 


r{uj TT'-f)/r{uJ tttttt) 


0.110 ±0.003 


0.0999 ±0.0026 ** 


r(0 TT'y)/r{(f) — > tttttt) 


(9.80 ±0.29) 10-3 


(8.3 ±0.6) 10-3 


r(pO ^ 7r7r)/r(pO ^ TT^) 


(2.52 ±0.03) 10^ 


(1.74 ± 0.27) 10^ 


T{uJ 7r7r)/r(co' 7r°7) 


0.169 ±0.009 


0.20 ±0.04 


r(pO ^ ^^)/r(pO ^ 7r°7) 


0.561 ±0.008 


(0.50 ±0.08) 


r{uj —>■ ri'-))/T{uj 7r°7) 


(6.62 ±0.17) 10-3 


(9.8 ±2.4) 10-3 ** 


r(0 ^ r/7)/r(0 ^ vr^) 


9.77 ±0.20 


10.9 ±0.3 ±0.7 



Table 6: Ratios of Widths. In the second data column, number written in plain style are directly 
extracted from the Review of Particle Properties, the boldface ones are derived from making 
the ratios of the accepted branching ratios |[T9l . The most significant results are flagged by **. 

As our fit procedure is global, it also allows to relate different decay modes of the same 
vector meson. These are provided in Table [6] with information extracted from the most recent 
RPP [19]. Taking into account the fit quality of the corresponding cross sections, the difference 
with RPP expectations should be considered significant. 

One should first remark the ratio r(p° — > 7r7) /r(p° — tttttt) which is found ~ 3.5 larger 
than expected from the ratio of branching ratios as given in [T9l, with a significance of ~ Aa. 
Less impressive but, nevertheless, significant is the corresponding ratio for the uj meson. 

More interesting, however, is the ratio V{uj r]'y)/T(uj 7r°7) which is found much 
more precise than (old) existing measurements and benefits from being derived by means of a 
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simultaneous fit to all recent relevant data sets. 



14 Conclusion 

The model developped in ^ succeeded in accounting for 17 decay modes of vector mesons 
and in providing a good simultaneous description of the e+e^ tt+tt^ cross section and of the 
dipion spectrum in r decay. The present paper has shown that it can be easily extended in order 
to provide also a successful description of the e+e~ 7r°7, e~^e~ 777 and e+e^ vt^tt+tt^ 
cross sections and account precisely for the properties of the r]/r]' 7r+7r~7 decays. Its 
application to e+e" KK decay mode is postponed to another study, because of the problem 
stated by d. 

This extension has been performed at the expense of only 2 more free parameters ci — C2 
and C3 which describe the amount by which the pure VMD assumption is violated. The precise 
numerical value of these are found close to the guess of Harada and Yamawaki [2J and amount 
to a small violation of the full VMD assumption in the Triangle Anomaly Sector and to a large 
violation for the Box Anomalies. On the other hand, we did not find any clear and unambiguous 
evidence for a high mass vector meson influence in e^e~ annihilations up to ~ 1.05 GeV. 

We have also shown that all experimental data, except for a data set on e^e^ tt^tt+tt" 
annihilation, are successfully described by a global fit within a unified (HLS) framework. This 
is clearly more constraining than only accounting for partial decay widths. Systematic uncer- 
tainties, especially correlated errors, have been carefully studied because of their importance in 
numerical estimation of some physics parameters. It has thus been shown that all data samples 
have already the expected absolute scale, except for the high statistics two-pion data sample 
from KLOE which has to be rescaled according to expectations. Then, all e+e^ data samples 
exhibit a good consistency with each other, including the KLOE sample which only marginally 
degrades the global fit quality. 

In summary the Extended HLS Model, supplied with the vector meson mixing mechanism 
defined in [T] provides a successful simultaneous description of all the data falling inside its 
scope. One should note that this model allows for the first simultaneous fit to all the low energy 
data not affected by scalar mesons or higher mass vector mesons. 

The last part of this study was devoted to results on partial widths of vector mesons as 
derived from our global fit. We have argued on the difficulty to firmly assess partial width 
values of vector mesons to some final state / relying only on e^e~ annihilation to this fi- 
nal state. We have preferred limiting ourselves to providing our estimations for products like 
V{y r+e~) X T{V /) or ratios of the form T{V fi)/T{V /2) where the uncertain- 
ties in estimating T(y e^e") cancel out. Taking into account the quality of our fits and the 
additional information provided by the physics underlying various processes, we believe our 
results are reliable and could supersede several older reported results. 
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Appendices 

A The Non-Anomalous HLS Lagrangian 

We remind flSl in this Appendix the HLS Lagrangian piece describing the photon sector 
(traditional VMD) : 



CvMD = te{l - .T,-dTx+ + i — {zA-'^- h)A ■K-dK+ + i—hA ■ d 



lag 



lag 



lag 

AZA 



-eagf^ 



Pi + 



V2 



Zv4>l 



A + lafy{5 + zy)A' + ^' 



(78) 



The parameter g is the traditional universal vector meson coupling constant ; the parameter 
a specific of the HLS model, is expected equal to 2 in standard VMD approaches but rather 
fitted to a ~ 2.3 2.5 lEIllMlEl. The parameter b in Eq. (TTSl) is 6 = a{zv — l)/6 where zy is 
the SU(3) breaking parameter of the Cy part of the HLS Lagrangian, while za = [/k/ fw]"^ = 
1.495 ± 0.031 (61 is the SU(3) breaking parameter of its Ca part Eia. 

A subscript / on the fields, standing for "ideal", affects the neutral vector meson fields. It 
indicates that the corresponding fields occuring in the Lagrangian are not the physical fields. 

On the other hand, the part of the HLS model involved in r decays writes [3J : 



a. 



d K- 



-p 



Za 



V2 



a r 



4 



za\Vus? + \Vud?] + -A^/^-za]\Vus?\w+ ■ W- + a^Vp+p 



4 



(79) 

plus the conjugate of the interaction term (the term, not displayed). This Lagrangian piece 
depends on the CKM matrix element = 0.97377 ± 0.00027 [[651, on g2 (fixed by its relation 
with the Fermi constant) : 

g2 = 2mw ^GfV2 , (80) 

on the universal coupling g and on the breaking parameters za and zy already defined. One 
should note, balancing the photon mass term in CyMD, a small mass term complementing the 
W mass of the Standard Model which could be removed by appropriate field redefinitions. At 
the T lepton mass scale one has [i65il : 

g2 = 0.629 (and e = 0.30286) . 
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B Effective Lagrangian Piece for e+e P7 

The VVP Lagrangian piece relevant for the P7 final state can be derived from the first 
Lagrangian in Eqs dH) ; it is given in terms of ideal fields by : 



' ^8 



(81) 

where the pseudoscalar fields are still the bare fields. 

As stated in the main text, the condition zwz^ = 1 is requested. Using this piece and Eq. 
(1781) . one can define an equivalent effective term which allows for a simpler derivation of the 
VP'^ couplings : 



(82) 



where now the pseudoscalar fields are the renormalized ones, Q is the quark charge matrix, A 
is the electromagnetic field and B' is given below. The effective piece of interest can further be 
written : 



H. 



a/3 



AVP — c u^j-i-yiiafj 




with B' 



vr ^ r]s ^ XT]o 
6 2VS 



vr ^ r]s ^ XT]o 
2 6^3 3^6 



(83) 



X 



3* - 71* 



in terms of renormalized ideal and pseudoscalar fields, renormalized after applying the SU(3)/U(3) 
breaking mechanism. One can now replace the ideal vector fields by their physical partners us- 
ing Eqs. (|22l) and gets : 



(84) 



where the sum extends over the physical (neutral) vector fields and the neutral pseudoscalar 
mesons. We have defined (A = \^/za) '■ 



i/^o = X 



1 + 3a 

3 + a - 2A/3 
6^3 
3 + a + A(3 



3^6 



= X 



3 — a 

1 - 3a - 2A7 
6^3 
1 - 3a + A7 



3^6 



HT 



/? + 37 
6 

7 + 3/3 + 2A 
6^3 
3/5 + 7 - A 



(85) 



X- 



3^6 
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where the dependence of the "angles" upon s has been omitted. 

The other Lagrangian piece contributing to e+e^ — ^ P7 annihilations is (see Eqs. dH) : 



Caap = B"e^-^^d^A,d^Ap K 

rO 



B" 



K 



vr ^ 5za -2 r]s ^ 5za + I x r]o 
6 3zA 6^3 ^za 3^6 



which provides the constant contributions to the amplitudes given by Eqs. (1411) . 
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C Integrated Cross Section For e+e 



Using the parametrization of Kuraev and Siligadze [43], the differential cross section is 
written : 



dx dy 



^^s^Gix, y)\Tsym + Tp\^ 



where : 



2 2 / 

G{x, y) = 4(x' - '^)(y^ - - [i - 2x - 2y + 2xy + 



r> 2 2 \ 2 

2m^ — rriQ \ 



(87) 



(88) 



The variables x and y are defined by (mo = m^o , m,r = ) by : 

s{2x + 2y-l) + ml 



The integrated cross section is : 



ale e — > tt tt tt 



■s+o ■ 

S-0 ■ 
0\ _ 



sil-2y) + ml 
s(l — 2x) + ml 



dx I dy — — 



with : 



Xr, 



i 1 



y min /max 



2(l-2x + xi 



mo(2m^ + mo) 



;i-x)(l -2a; + 



{2ml -ml). 



+/- 



ml 



(x' -)(l-2x + 

s 



mo(2m^ - mo) 



1 - 2x 



mo(2m^ + mo) 



1/2- 



(89) 



(90) 



(91) 
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Figure 1: Annihilation process e+e n^^. The errors plotted combine the reported system- 
atic and statistical errors in quadrature. The insets magnify the top part of the to region on the 
one hand and the (f) region on the other hand. "CMD-2 2005" refers to 1361, "SND 2003" to 
Il39l and "SND 2000" to L40J. 
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Figure 2: Annihilation process e+e" 777. The errors plotted combine the reported systematic 
and statistical errors in quadrature. The insets magnify the cu region and the (p region. "CMD-2 
2005" refers to ^ (with r] ^ 27), . "CMD-2 2001" to [66] (with r] 37r°), "CMD-2 1999" 
to [37] (with r] Tc+n-ir^), "SND 2006-7a/b" refer to the two data sets published in JMl 
(with 7] 37r° and r] tt+tt^tt") and "SND 2000" to the data set published in [|40l (with 

V ^ 27). 
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Figure 3: ISR data [13] on the e+e^ vr+vr" cross section collected by the KLOE Collab- 
oration at DA$NE with the fit superimposed. The plotted errors are the square roots of the 
diagonal terms in the full error covariance matrix (see text). The inset magnifies the p — to 
region. Downmost plot shows the residual spectrum (with same units) coming out from the 
global fit performed with all Novosibirsk e^e^ tt+tt^ and e+e^ {t^^ Iv)! data. 
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Figure 4: Plots of CMD-2 data (left figures) and SND data (right figures) superimposed with 
the cross section predicted from running the global fit with the ND and CMD data as single 
3-pion data sets. Top figures show the case in the u peak region, downmost figures show the 
peak region. See text for reference and comments. 
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Figure 5: Fit performed using CMD-2, ND and the old CMD data (see text). Top figures show 
the uj and peak regions, the downmost plot shows, magnified, the region outside the peaks. 
The measurements of the former DM 1 data set have not been used in the fit, but are nevertheless 
displayed in this plot. See text for references and comments. 
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Figure 6: The isospin breaking parameter a{s). 
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Figure 7: The isospin breaking parameter P{s). 
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Figure 9: Dipion invariant mass spectra in the r/' 7r+7r^7 decays from Argus fSG'l, Lepton F 
[|58l and Crystal Barrel [59J experiments. The last plot shows the photon momentum spectrum 
in the decay r] -k^-k^^ from [|53l . The curves superimposed are predictions, not fits. 



58 



7 




59 



